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Abstract 

Some invariant sets may attract a nearby set of initial conditions but nonetheless 
repel a complementary nearby set of initial conditions. For a given invariant set X C W 1 
with a basin of attraction N, we define a stability index &{x) of a point x € X that 
characterizes the local extent of the basin. Let B e denote a ball of radius e about x. 
If cr(x) > 0, then the measure of B e \ N relative the measure of the ball is 0(e' CT ^'-"), 
while if (j{x) < 0, then the measure of B e n N relative the measure of the ball is of this 
order. We show that this index is constant along trajectories, and we relate this orbit 
invariant to other notions of stability such as Milnor attraction, essential asymptotic 
stability and asymptotic stability relative to a positive measure set. We adapt the 
definition to local basins of attraction (i.e. where N is defined as the set of initial 
conditions that are in the basin and whose trajectories remain local to X). 

This stability index is particularly useful for discussing the stability of robust hete- 
roclinic cycles, where several authors have studied the appearance of cusps of instability 
near cycles that are Milnor attractors. We study simple (robust heteroclinic) cycles 
in M 4 and show that the local stability indices (and hence local stability properties) 
can be calculated in terms of the eigenvalues of the linearization of the vector field at 
steady states on the cycle. In doing this, we extend previous results of Krupa and Mel- 
bourne (1995,2004) and give criteria for simple heteroclinic cycles in M. to be Milnor 
attractors. 
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1 Introduction 

For many choices of smooth vector field / : IR n — > IR n , the system on x £ W 1 

x = f(x) (1) 

has a small subset of M n (an attractor) that attracts a large set of initial conditions; these 
attractors are important for understanding the long term behaviour of trajectories of the 
system. In this paper we explore the local attraction structure of invariant sets, while in 
the latter part we focus on a particular class of examples - attracting heteroclinic cycles. 
More precisely, an invariant set is asymptotically stable if it attracts all nearby points; many 
systems are found to possess invariant sets that are not asymptotically stable, but that are 
attractors in a weaker sense (e.g. in the sense of Milnor |2U]). 

Now consider £1, . . . , £ m to be hyperbolic equilibria of ([I]). A set of connecting trajectories 
W u {^j) PI W s (£j+i) ^ 0, j — 1, . . . , m, £ m +i = £i, is called a heteroclinic cycle between these 
equilibria. It has been shown that heteroclinic cycles can be robust (persistent to small per- 
turbations) if / is constrained to be symmetric with respect to certain group representations 
[T5| HI [22] , or if / is constrained to preserve certain invariant subspaces [32] . Heteroclinic 
cycles that are not asymptotically stable may often be observed to be apparently stable in 
computations. To explain this, weaker notions of stability for heteroclinic cycles sets were 
introduced in [191 13 EES] - they do not require attraction in a full neighbourhood of the in- 
variant set; they may even be repelling in a region that is typically cusp-shaped in Poincare 
sections to the cycle. The papers [TjJJ |TT] define a heteroclinic cycle to be essentially asymp- 
totically stable (e.a.s.) if it attracts almost all nearby trajectories, and they define it to be 
almost completely unstable (a.c.u.), if it attracts almost no nearby trajectories. However, as 
shown in [13] these definitions are not mutually exclusive. (Brannath [9] similarly discusses 
e.a.s. using the notion of relative asymptotic stability from Ura |23j.) 

The paper is organized as follows: in Section |2] we discuss various definitions of stability, 
and we relate them to the notion of Milnor attractor and the local geometry of the basin of 
attraction. We introduce a stability index that characterizes the local geometry of the basin 
of attraction. After proving some basic properties about this invariant of the dynamics, 
we generalize to a local stability index that is the limit of stability indices of local basins of 
attraction. In Section [3] we discuss the structure of heteroclinic cycles and describe the geom- 
etry of local basins of attraction by way of the local stability index and (Poincare) surfaces 
of section. We show, under certain assumptions, that the stability index of a connecting 
trajectory is the stability index on a surface of section. 

Section H] computes the stability indices for robust heteroclinic cycles in M 4 ; we employ 
the classification of simple cycles in M 4 into Types A-C by Krupa and Melbourne in a 
series of papers [Ml H3 [TS] and calculate the stability indices of the connections in terms 
of eigenvalues of the linearization at equilibria in the cycle. Finally we discuss some of the 
limitations and possible further uses of stability indices and related concepts in Section [5j 
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2 Attractors and the stability index 



Various definitions of attraction of invariant sets have been introduced [9l [T7J [191 122] to 
describe sets that are not asymptotically stable but that are nevertheless attracting in some 
sense. We review these notions and relate them to Milnor's notion of a measure attractor 

2.1 Notions of attraction for invariant sets 

In this section we consider a smooth flow *&t( x ) on Two very general notions of attraction 
are the Milnor and weak attractors discussed in [20] and [6] respectively. For an invariant 
set X C M. n we define the (global) 
basin of attraction of X to be 

B(X) = {i£ R n : u(x) C X} 

where u(x) = Qrx) {^O^) '■ t > T} is the w-limit of x. The following defines attraction 
properties of X in terms of this basin. We use £(■) to denote Lebesgue measure on M n . 

Definition 1 |2J/ We say a compact invariant set X is a weak attractor if £(B(X)) > 0. 
We say a compact invariant set X is a Milnor attractor if it is a weak attractor such that 
for any proper subset Y C X that is compact and invariant we have 

£(B(X)\B(Y)) > 0. 

We do not assume transitivity of X (a dense orbit); indeed, the main examples we will 
consider later on are heteroclinic cycles that are not transitive. Note that any weak attractor 
contains a Milnor attractor [6j Lemma 3.2]. There are various examples of robust heteroclinic 
cycles (e.g. [T9J [9J [13]) that are Milnor attractors, even though they are not asymptotically 
stable. Let d(-, •) denote the Hausdorff distance between two sets, let 

B e (X) := {x G W n : d(x,X) < e} 

denote the e-parallel body of X, and let D c denote the complement of D in M. n . 

Definition 2 J7P) / We say a compact invariant set X is essentially asymptotically stable 
(e.a.s.), if there is a set D such that for any open neighbourhood U of X and any e > there 
exists an open neighbourhood V C U of X such that: 

(a) Ifx eVnD c then $<(x ) G U for all t > and lim^oo d($ t (x ),X) = 0, 

(b) £(V fl D C )/£(V) > 1 - e. 
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Intuitively, if X is e.a.s. one might expect that it attracts "almost all" nearby trajectories, 
while [T7j says X is almost completely unstabl^j] if it attracts "almost none" of them. How- 
ever, these definitions do not formalise these intuitive categories very well; as highlighted in 
[T3] , they are not mutually exclusive and so may yield classifications that are not intuitively 
helpful. Another useful definition is that of [23] which is used in [9]: for this we consider a 
set iV C r. 

Definition 3 JMBj We say a compact invariant set X with X G N is asymptotically stable, 
relative to (a.s.r.t.) N if for every neighbourhood U of X there is a neighbourhood V of X 
such that for all initial x G V fl iV we have <&t{x) G U for t > 0, and u(x) C X . 

In fact, Brannath [9] interestingly suggests that the authors of [191 HZ] had the following 
definition in mind for e.a.s., but we name it differently to distinguish from the original 
definition in [T9l. 



Definition 4 (Adapted from f^) We say a compact invariant setX is predominantly asymp- 
totically stable (p.a.s.) if there is an N such that X is asymptotically stable relative to N 
and 

emx)nN) 

e->o £(B e (X)) 

We now give a result that relates these concepts of attraction. 

Theorem 2.1 Suppose that X is a compact invariant set for a continuous flow 

(a) X is p.a.s. =>• X is e.a.s. 

(b) X is e.a.s. =>■ X contains a Milnor attractor. 

(c) X is e.a.s. there is an N with £(N fl A) > for any neighbourhood A of X , such 
that X is a.s.r.t. N . 

Before proving this theorem, we give a useful lemma that will be used in the proof. For 
any measurable set N we define the density of N at x to be 



. , ..Ax) n N) 
Fix) = lim v n ) K J t - s 1 

and recall that the Lebesgue Density Theorem [TT] states that for ^-almost all x G N we 
have F(x) = 1. In such a case we say that x is a point of Lebesgue density for N. 



1 A flow- invariant set X is called almost completely unstable (a.c.u), if there is a set D and an open 
neighbourhood U of X such that for some e > there exists an open neighbourhood V of X, V C U, such 
that (a) for x £ V \ D there exists a t > with $ t (x ) £ U; and (b) £(V n D C )/£(V) > 1 - e. 
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Lemma 2.1 Suppose that N has positive measure and Y be any closed and bounded subset 
of N with zero measure. Then for any e > one can find an open set V containing Y with 

£{Vf]N)/£(V) > 1 - e. 

Proof: Although Y need not contain any points of Lebesgue density for N, there is at least 
one point x & N of Lebesgue density, and so we choose 5 > such that 

£(B s (x)nN)/£{B s {x))>l-^. 

Now let V = Bs(x) U B V (Y). Because of outer regularity of £, r\ can be chosen small enough 
to ensure that £{B rj {Y)) is as small as desired, and hence the result holds. QED 

With a slight modification of the argument, one can assume that V is connected and 
open in the statement of the above Lemma; however it may be very far from being a ball in 
terms of the relationship between diameter and volume of the set. 

Proof: [of Theorem 12. lj For (a) suppose that X is p.a.s. and let N be a set for which X is 
a.s.r.t.. By Definition HJ for any e > there exists a 5 > such that 

£(B s (X)nN) 
£(B S (X)) 

for all 5 < Sq. In Definition [2] we set N = D c and V = B$(X) (where 5 is sufficiently small 
so that Bs(X) C U), and we prove that X is e.a.s.. For (b), note that this follows because 
iV = D c is a subset of the basin of attraction of X and has positive measure as £(V fl N) > 
for some set V. Hence it is a weak attractor, and contains a Milnor attractor [6]. Finally, 
for case (c) suppose firstly that X is e.a.s., then it is stable relative to the set N = D c and 
£(N n A) > for any neighbourhood A of X. The converse for (c) follows similarly, on 
applying Lemma 12.11 QED 

There are examples that show that, in general, converses of (a,b) do not hold; for a 
counterexample to the converse of (a) we refer to [13] who present heteroclinic cycles that are, 
in our terminology, e.a.s. but not p.a.s.. For a counterexample to the converse of (b), there 
are "unstable attractors" [7], though only for a weaker assumption - that $t is a semiflow. 
These "unstable attractors" are Milnor attractors that have zero basin measures within a 
small enough neighbourhood of the attractor. It is not clear whether the converse of (b) is 
true for flows (possibly subject to some smoothness assumptions). Note that Theorem 12.11 
is a generalization of comments already made in [9J pl369] which assume iV to be an open 
set. There are examples of heteroclinic cycles that are not asymptotically stable relative to 
any open set, but that do seem to be asymptotically stable relative to a positive measure 
"riddled" set 0. 
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G(X)>0 



G(X)<0 



Figure 1: Schematic diagram illustrating how the stability index o~(x) of a point x G X 
relates to the local geometry of the basin of attraction of X (shaded region). For cr(x) > 0, 
the measure of points in a ball of radius r that are in the complement of the basin goes to 
zero, relative the measure of the ball, as r' "^'. For o~(x) < 0, this estimate applies to the 
basin itself. 



2.2 Geometry of global basins: the stability index 



We suggest that it is useful to distinguish between different local geometries for e.a.s. sets. 
To this end, consider X an invariant set in IR n and let N = B(X) denote its (global) basin 
of attraction; we assume that the flow $ t is smooth. Pick a point define 



£ e (x) 



(x) n N) 



£(B £ (x)) 
and note that < S e (x) < 1. 

Definition 5 For a point x G X we define the stability index of X at x to be 



(2) 



which exists when the following converge 

"ln(E £ (x)) 



<J-[X) 



lim 

e-S-0 



ln(e) 



<7_|_(XJ _ C- 



a + (x) 



x 



lim 

<E->0 



ln(l - E e (x)) 
ln(e) 



We use the convention that cr_(x) = oo if there is an €q > such that £ e (x) = for all 
e < e 0; and o~ + (x) = oo if E e (x) = 1 for all e < e . Note that <r±(x) > and so we can 
assume that a(x) G [— oo, oo]. 

The stability index may not exist at certain points in X (an example is given in Section |5]), 
and may vary throughout X when it does exist. Figure [1] illustrates how the local geometry 
of the basin relates to the sign of a(x) for a point x G X. Note that o~(x) = +oo is 
the "strongest" form of local stability while — oo is the "weakest". The following Lemma 
characterizes some basic properties of the index: 
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Lemma 2.2 Suppose that a(x) is defined for some x G X C M. n ; then the following hold: 

(a) If one ofa±(x) converges to a positive value then the other converges to zero (i.e. only 
one of a + (x) and cr_(x) can be non-zero). 

(b) If a(x) = c > then 1 — S £ (x) = 0(e c ) (and in particular T, t (x) — > 1 as e — > 0). 

(c) If a(x) = — c < then S e (x) = 0(e c ) (and m particular Tl e (x) — > as e — )■ Oj. 

Proof: For (a) note that if cr_(x) > then lim e ^ S e (x) = 0; this implies that 1 — S e 
converges to 1 as e — >■ and so cr+(x) = 0. The other case is argued in a similar way. (b) 
follows on noting by (a) that c = a(x) = cr+(x) > and cr-(x) = 0. Hence we have from the 
definition of a+(x) that 1 — S e (x) = 0(e c ). A similar argument gives (c). QED 

The main result in this section is the following; this can be generalized to cases where 
a(x) is measured relative to any measurable invariant set N. 

Theorem 2.2 Suppose that N = B(X) is the basin of a compact invariant set X for a C 1 - 
smooth flow 3>t(x). Then for any x the index a(x) is constant on trajectories, whenever it is 
defined. 

Proof: Fix x G X such that o~(x) is defined, and pick any t > 0. Let (f)(x) = 3>t(x). Because 
is a C 1 diffeomorphism, one can find an i] > such that there is an L > 1 and an M > 1 
with 

-j- < det(D(f>(y)) < L, -j^ < \\Dtp(y)\\ < M (3) 

for all y G B v (x), where D denotes the derivative (Jacobian) of the map. We assume that 
L, M are chosen so that the same inequalities are satisfied by Dcf)" 1 for z G B v (<f)(x)). As a 
consequence of this, one can find rj' with < rj' < r\ such that 

B £/M (x) C <jr\B € {<i>{x))) C B Me {x) (4) 

for any e < rj'. Writing Xn{v) as the indicator function for and y G B v (x) we have 

A = £( J B e (0(x))niV) = / XN(y)d£(y) 

J yeBt(<f>(x)) 

= [ XN{z)det{Dcj ) - 1 {z))di{z) 
Jzet-HBtMz))) 

where in the last line we have substituted y = <f>(z) and we have used the fact that is 
invariant, so that 

X N {y) = 1 & Xn{<P{v)) = 1. (5) 
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Hence, using (|3 

1 

L 



XN{z)di{z) < 



z&B e/M {x) 



< 



Xn{z) d£(z) 

z&<t>-i(B e {4>{w)) 

X n{z) &et(D<p-\z)) 



A 



< L 

< L 



Xn{z) d£{z) 



*e^- 1 (B«(*(*))) 

Xn(z) d£(z) 
zeB eM (x) 



meaning that for all e < rf we have 

j£(B e/M (x) HN)< £{B e {4>{x)) n iV) < L£(B eM (x) n iV). (6) 

This means that from (j2j), there is a A" = LM" such that for all small enough e we have 

—Y, e/M (x) < E e (0(ar)) < i^S eM (x) 

(we have used the property that £{Bme{x)) = M n £(B e (x)) for any e > and x). Hence we 
have 



In e - In M 
In e 



ln(S e/M (x)) In AT ln(S €/M (x)) In AT _ In (±X e/M (x)) 
ln(e/M) lne 



In e 



In e 



< 



< 



In e 
lne + InM 

In e 



and taking the limits as e — > we have 

cr_(x) < cr_(0(a;)) < (T-(x). 



In e 

ln(S e (0(x))) 
In e 

ln(A-S eM (x)) ln(S eM (a;)) + lnK 



In e 



In e 



ln(S eM (x)) + lnA- 



ln(eM) 



In e 



(7) 



A similar argument on substituting by its complement gives cr + (x) = a + (4>(x)) and hence 
the value of a(x) is constant along trajectories of $ t . QED 

Note that this argument works for any C 1 -diffeomorphism <fr for which A^ is invariant, 
meaning the result can be used to show that cr(x) is invariant under C 1 -conjugation - it is an 
invariant of the dynamics. Note also that although cr(x) is constant on a given trajectory, it 
may depend on which trajectory is chosen. 

The stability index can be used to determine e.a.s. and p.a.s. by the following theorem. 
However, converses of the following theorem are not expected to be true in general as c{x) 
may be negative on a "lower dimensional" set of trajectories within X, or may not converge. 
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Theorem 2.3 Suppose that for all x G X the stability index o~(x) G [—00, 00] is defined. 

• If there is a point x G X such that —00 < a(x) then X is essentially asymptotically 
stable (e.a.s.), and contains a Milnor attractor. 

• If there is a c > such that c < a(x) for all x G X then X is predominantly asymp- 
totically stable (p.a.s.). 

Proof: (a) The fact that —00 < cr(x) implies in particular that N = B(X) contains a set of 
positive measure, and so by Theorem 12.1( c) it is e.a.s.. By the Theorem 12.1( b). X contains 
a Milnor attractor. (b) Note that the basin of attraction N of X is such that for any 5, 
S e (x) > 1 — 5 for all x, and some e depending on x. By compactness of X one can choose 
an e small enough that i(B e (X) n N) > (1 - 6)£(B e (X)), implying p.a.s. of X. QED 



2.3 The local stability index 

While Definition considers the global basin of attraction, the stability index a(x) can be 
adapted to provide a useful concept from purely local properties of the attractor. We define 
the 5-local basin of attraction to be the basin of attraction of X relative to Bg(X); that is, 

B S (X) := {x : u(x) C X and $ t (ar) G B S (X) for all t > 0}. (8) 

Note that B&{X) is forwards, but not necessarily backwards invariant under the flow. The 
limit of the stability index for points relative to the 5-local basin as 5 — > is called the local 
stability index o~i oc (x) for X. More precisely, we define 



£(B e (x)nB s (X)) 
(*)) 



S f Ax) 



and for a point i6lwe define the local stability index of X at x to be 



(9) 



which exists when the following converge 



cr loc (x) := limlim 

8^0 e^O 



ln(S e ^(o;)) 
ln(e) 



dine +(x) := lim lim 



ln(l - S £ ^(x)) 
ln(e) 



with the same conventions as before. The definition works for discrete time systems as 
well as continuous time, without further modification. Note that the local stability index is 
computed for 5 small and fixed before taking the limit as 5 — > 0. 

One can weaken the assumptions in Theorem 12.21 to give the same conclusion; the critical 
step is that if r\ < S and we assume that G B n (X) for < s < t then (JSJ) still holds- 

and by continuity of one can choose a small enough e that x and <f>(x) are guaranteed to 
be in such a B V (X). 
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2.4 Stability indices for sections to the flow 

Suppose that $t(x) has an attractor X C M. n and pick a point x £ X. Let S C 1™ be a 
smooth n — 1-dimensional subspace containing x that is transverse to the flow at x. One can 
relate the stability index cr(x) or a\ oc (x) to the stability index for the dynamics defined by 
the return map F on S as follows. 

Theorem 2.4 Suppose that X is invariant for a ^-smooth flow $t(x) and that N is a 
(local) basin for X . Suppose that S is a codimension one surface that is transverse to the 
flow at x; then a(x) can be computed relative to the intersection of N with S on substituting 
E e (x) by 

e s (B e (x)nNns) 
^' s[x) i s (B e (x)nS) ■ 

Proof: Let N = B(X); the argument for local basins will be similar. Note that iV is 
invariant implies that it is a union of trajectories. We consider local coordinates in M. n 
near x that are the coordinates in S and time. Pick any small e > 0; by simple geometric 
arguments (i.e. you can always put a cylinder in a larger sphere, and a sphere in a larger 
cylinder) there is a constant K > 1 such that 

B e/K (x) C {B e (x) HS)x [-e,e] C B eK (x). 

Using the product structure of Lebesgue measure we have 

£{B e/K {x)) < 2e x £ s {B e {x) n S) < £(B eK (x)) 

with a similar inequality for £s{B e {x) H N H S). Hence 

• e/K {x) n N) < £ s (B e (x) DNDS) < £(B eK (x) n N) 



meaning that 



[B eK (x)) ^ £ s (B e (x)nS) £(B £/K (x)) 



£ s (B e (x) n 5) 

and as before (for the flow) S e (x) satisfies the inequalities 

1 



Z e/K (x) < £ e , 5 (x) < (2K) n E eK (x) 



(2K) n ' 

where we have used the fact that £(B t x(x)) = (2K) n £(B e / K (x)). In particular, the scalings 
of these quantities are the same as e — > 0. QED 

Theorem 12.41 implies, for example, that if there is a return map for the flow on S then the 
stability index of trajectories for a flow can be computed by examining the stability index 
for the intersection of the basin with a suitable surface of section. 
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3 Robust heteroclinic cycles 



Suppose that T is a finite group acting orthogonally on M n , and that / : M. n — > M. n is a 
T-equivariant vector field, i.e. 

f( JX ) = 1 f(x), for all 7 G T. 
Let j — 1, . . . , m, be hyperbolic equilibria for 

i = /(a;) 

with stable and unstable manifolds W s (^j) and respectively, and let Sj = W u (£j) fl 

W s (£j+i) 7^ be connections between ^ and £7+1, where £ m +i = £1; then the group orbit X 
of the equilibria and the connections 

X = clos ({7Sj : j = 1, . . . , m, 7 G T}) 

is called a heteroclinic cycle. Recall that for a group T acting on W 1 , the isotropy of the 
point x G W 1 is the subgroup 

E x = {7 G T : 7X = x} 

while for a subgroup S C T, a fixed-point subspace of E is the linear subspace 

Fix(E) = {1 6 K" : (jx = i for all a G E}. 

In the absence of symmetry or other constraints, a vector field with a heteroclinic cycle 
is structurally unstable, i.e. there are arbitrarily small perturbations of / to g, such that 
the heteroclinic cycle does not exist for the vector field g. For symmetric vector fields, 
heteroclinic cycles may be robust, as long as each connection is robust within some invariant 
subspace and only symmetric perturbations are allowed [T5] . 

3.1 Local structure: eigenspaces and simple cycles 

A sufficient condition for a cycle X to be structurally stable (or robust), is that for all j there 
exists a subspace Pj such that Pj = Fix(Ej) for some Ej C T, Sj C Pj, is a sink in 
Pj. Denote Lj = Pj fl Pj-i- We denote the isotropy subgroup of points in Lj \ {0} by Tj. 
Note that X QY, where Y is a linear subspaces of the inner product space X, denotes the 
orthogonal complement to Y in X. 

If X is a structurally stable heteroclinic cycle then the eigenvalues of (df)^. can be divided 
into four classes: 

• Eigenvalues with associated eigenvectors in Lj are called radial, the maximal real part 
of radial eigenvalues being —rj. 

• Eigenvalues with associated eigenvectors in Pj-% Lj are called contracting, the max- 
imal real part of contracting eigenvalues being —Cj. 
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• Eigenvalues with associated eigenvectors in Pj Lj are called expanding, the maximal 
real part of expanding eigenvalues being ej. 

• The remaining eigenvalues are called transverse, the maximal real part of transverse 
eigenvalues being tj. 

The heteroclinic cycle X e R 4 \ {0} is called a simple robust heteroclinic cycle (in M 4 ) 
PJ if for all j: 

• All eigenvalues of (df)^ are distinct, dimPj = 2 and X intersects with each connected 
component of Lj \ {0} in at most one point. 

For simple cycles, Lj and the three remaining subspaces are one- dimensional, hence 
there is a unique real eigenvalue of each type. Moreover, for simple cycles either Tj = 7% 
and T,j = Z 2 for all j or Tj = J\ and E 3 - = Z| for all j (see Proposition 3.1 in [18]), and each 
simple cycle is of one of three types discussed by [TB"] . 

Definition 6 Suppose that X is a simple heteroclinic cycle that is robust for a vector field 
in R 4 with a finite symmetry group. We say 

• X is of Type A if Sj = Z 2 for all j . 

• X is of Type B if there is a subspace Q of M 4 with dim(Q) = 3 such that Q = fix(S) 
for some S CT and X C Q. 

• X is of Type C if it is neither of Type A nor of Type B. 

The work of [18] goes on to differentiate between four varieties of Type B cycles (denoted 
by B±, B£, B± and B^) and three varieties of Type C cycles (denoted by C-f , and C±), 
depending on the number of equilibria involved in the cycle and action of the group V. In 
Section H] we examine the stability of cycles in M 4 using Poincare maps, where the structure 
of the maps depend on the type of cycle. 



3.2 Local stability for heteroclinic cycles 

For a heteroclinic cycle X comprised of one- dimensional connections Sj,j = l,...,m(sj is 
the connection from to £_,•), its local attraction properties are described by the set of 
stability indices of the trajectories 

a = (ax, . . . , a m ). 

where 

Oj = a{x) 

for x an arbitrary point on Sj. The following lemma will be useful later on: 
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Lemma 3.1 Let a simple heteroclinic cycle X be comprised of one- dimensional connections 
and suppose that — oo < Oj for some j . Then X is a Milnor attractor. 

Proof: This follows from the fact that — oo < o~j(x) implies that £(B(X)) > and so X is 
a weak attractor. Since no invariant subset of the cycle can be a Milnor attractor, X must 
itself be a Milnor attractor. QED 

In what follows, we calculate local stability indices for different types of simple robust 
heteroclinic cycles in K 4 . Following [131 HE], to examine stability we construct a Poincare 
map in the vicinity of the cycle. 



3.3 Stability indices for return maps 

Section [3] gave definitions for radial, contracting, expanding and transverse eigenvalues of 
the linearization (df)^.. Simple cycles in R 4 will possess a single eigenvalue of each type. Let 
(u, v, w, z) be local coordinates near £j in the basis of the four associated eigenvectors, and 
Bg(£j) be a neighbourhood of £j defined as 

= {(u,v,w,z) | max(|tt|, \v\, \w\, \z\) < 5} 

where 5 is small, denote by (u,v,w,z) the scaled coordinates (u,v,w,z) = (u,v,w,z)/5. In 
the (u, v) plane we will also employ plane polar coordinates (r, u = r cos 6 and v = r sin 6. 
For a small 5, the vector field / can be linearly approximated in B$(£j). We assume that 5 
in OH]) is sufficiently small, so that 

B S (Q C B S (Q for all j. 
In Bg(£j) we consider the linearised system (CQ) 

ii = —rju 
w = ejW 

Z = tjZ 

This gives an accurate approximation of the nonlinear flow as long as the linear system has 
no low order resonances. 

The connection Sj is tangent to the subspace u = v = z = 0. The heteroclinic connection 
to £j lies in Pj-i where local coordinates are u and v. For 

Hj OUt ^ = {(u, v , w, z) : |i>|, \z\ < 1, w = 1}, 
Hj = {(r, 9, w, z) : r = 1, |iw|, \z\ < 1}. 
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the first return map <pj : 

H (in) ^ H (out) ig defined 

near each equilibrium in = Qj PI 

Hf n \ where Qf n) = {(z ,w)\ \z\ < |ty|* 3 ^ ej ')|j. For each connection we define connecting 



diffeomorphisms ipj : Hj OUl> — > Hj™[ and their compositions gj = ipj o : ifj m - ) — >• Hj 1 ^. 
The Poincare map is the composition g = g m o . . . o g x : — >• if £ m ' . 
As shown in [T3l HE] , at leading order the maps have the formal 

4>j(u,v,w, z) = (u\w\ rj/e: > ,vw C:i/ej , l,z\w\~ tj/ej ) (11) 

and 

ipjiu, v, w, z) = (1, 9 + (3jU, a^v + aj 2 z, oij z v + (Xj 4 z). (12) 

For these maps only (w, z) coordinates are important [131 [16], an d restricting to these coor- 
dinates the map gj : R 2 — > R 2 we have 

gjiw, z) = (a^w ^^ + a j2 z\w\~ tj/ej ,a j3 w Cj/e: > + a j4 z\w\~ t J /ej ). (13) 

For cycles of Type A, generically ^ for all j = 1, . . . , m, k = 1, 2, 3, 4. For cycles of 
Type B, atj 2 = aj 3 = and cx^a^ ^ 0. For cycles of Type C, aj l = aj 4 = and oij 2 oij z ^ 0. 
Thus, for a point x = X R h[ on the cycle X we have associated the map 

9 = 9m° g m -\ o . . . o g x : R 2 — > R 2 . (14) 

We will also use the notation gi^ = g\ o o . . . o g x o g k . 

Similarly to the local stability index for a point x G X (see Section 12.31) . we define a 
stability index a 9 for the map (fT4"j) . Note that -B e is the ball of radius e in R n centered at 0, 
and we define 

B 9 S := {x : xe R™, |# iJfe x| < 5 for all < I < m - 1, k > 0}, (15) 

to be the 5-local basin of attraction of in R n for the map g (T14T) . The local stability index 
is defined to be 



a 9 := at - a 9 



where 



o 9 _ := limlim 

<5-s>0 e^O 



ln(£ e 



a 9 , := lim lim 



ln(l - S e 
ln(6) 



2 Strictly speaking, the maps are defined for \z\ < K(l ~ <5)|w|* J / ej , where K is a constant and 5 is small 
[13) . For simplicity, we ignore K and S, because for small z and w they do not enter into asymptotically 
significant terms. 

3 The maps for negative w are defined as follows. Consider a neighbourhood of a point £j. Two heteroclinic 
connections enter this neighbourhood, — > £j and 7 < - m - > ^-i — ► £j, where 7 ( m ) is any symmetry, satisfying 
£ Xj and 'y( m ) ^j-i- Two heteroclinic connections exit the neighbourhood: £j — > and £j — ► 
7 (ou*)^. +lj w here 7 ( ou *) is a symmetry, satisfying ^ out ^ G Tj and 7^"^ £ Ej. The local map ^ : i?] m) -> 
jj(out) j g (iefixieci only for u; and v of particular signs, say w > and u > 0. For w < 0, the local map acts to 
Hj° ut \ where ffj ou *' — ^(o«t)^( otlt ). £ or ^ < 0, it is defined in j( m ) Hj ln \ Due to existence of the symmetries 
ry{in) an( ^ ry(out) ^ we c&n consider <pj for arbitrary w and w: by applying these symmetries the local map can 
be defined for w and v of arbitrary signs. 



14 



with 

^ _£(B e (o)nBI) (16) 



,(o)) ■ 

Because of the asymptotic independence of the return map on two of the coordinates, we 
can effectively reduce the computation of the stability index for heteroclinic cycles in R 4 to 
a calculation on a section to the cycle in R 2 . 

Theorem 3.1 Let g : R 2 — > R 2 [Lfy be the map associated with a point x = X fl H[ m \ 
where X is a simple heteroclinic cycle in R 4 . Then 

^locO) = o- 9 . 

Proof: First, we prove that o"i oc (x) > a 9 . Denote Q = [— e, e] 2 x B 9 . If y = (u,v,w,z) G 
B t (x) \ Q, then there exist j and k, such that \gj,k{y)\ > 5. For small 5, a trajectory near the 
heteroclinic cycle is approximated by the maps <pj (TTTT) and ipj (Tl2l) . where the coordinates 
(w, z) are independent of (u,v). Hence, for the point (u, v, w, z), which is the {k + l)st 
intersection of §t{y) with Hf n \ we have \(w, z)\ > 5/2. Hence y ^ i3<w 2 (X), and therefore 



6/2\ 



(x)n^ (5/2 (X))<4e 2 £( J B e ni3 9 



implying that o"i oc (x) > a 9 . 

Second, we prove that o"i oc (a;) < a 9 . If y = (u,v,w, z) G £> £ (.t) fl Q, then for any 
intersection (l,9i,w,z) (here we are using polar coordinates in the (u, v ) plane, 9i is the 
difference between 9 and #0, the distance of the point from the cycle is denoted \(9i,w,z)\) 
of $t{y) with Hj tn \ \(w, z)\ < 25 holds true. Together with ffTTj) and ffl2l) . it implies that in 
any intersection 9\ < (35 r , where /3 = max (max.,- (/?,,' ), 1) and r = min(mm j(rj/ej), 1). Thus, 
we have proved that for any intersection y of with Hj for any j we have 

d(y,X) = \(9 1 ,w,z)\<2p6 r . (17) 

If a trajectory is close to the heteroclinic cycle, then there exist constants Kj such that 
d(sj, $t(x)) < Kjd(sj, $T J fc (^)) in the interval < i < (Here 2}^ is the time when 

$t(x) crosses Hj and Tj_i^ is the time when it crosses iij™ .) Denote = /^max^i^). 
In the vicinity of £j we can consider the system (TTTJT) using the approximated map /, hence 
(|T7|) is satisfied at the points of intersection, which implies that 



d(® t (y),X) < 2K5 r for all t > 0. 
Hence y G £>2if<5'' (^0 , and therefore 

£{B e (x) n 2WP0) > e 2 £(5 e n £f ). 

Since K > and r > do not depend on 5, this implies cri oc (x) < a 9 and therefore 
<Tioc(z) = a 9 . QED 
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4 Stability indices for heteroclinic cycles in M 4 



In this section the stability indices for the connections of simple robust heteroclinic cycles 
in R 4 are calculated in terms of ratios of eigenvalues aj = Cj/ej and bj = —tj/ej, 1 < j < m. 
We do this relative to the classification of simple heteroclinic cycles in M 4 of Definition 
and [TS]. Here only statements of the main theorems and a sketch of the proof of the main 
theorem for type A cycles are presented. The complete proofs are given in the Appendices. 

Using the maps {gj , j — 1, . . . , m} from the previous section, we define return maps on 
sections to any of the connections 

g(j) . fff») 
3 J 

to be g& = gf,-_i ° . . . o gi o g m o . . . o gj, let 



and 



a 3,± = (T ± 

where a± are defined in Section 13.31 Note that in the previous section we used 

9 = 9 {1) 

and that aj will effectively give the stability index for the connection that intersects H ■ . 



4.1 Type A cycles 

Consider the map g = g m o . . . o g x ; ]R 2 — y ]R 2 where 

gj(w,z) = {gj{w,z),g z j {w,z)) = {A jW a ^ + Bj\w\ bj z, CjW " 3 + Dj\w\ h 'z). (18) 

For generic cycles of Type A we can assume AjBjCjDj ^ and AjDj ^ BjCj for all j. 
Recall that — gi o . . . o g x o g k . 
Let us define 

Q(r,e) = {{w,z) : \(w,z)\<e, \z\ > \w\ l+r , \w\ > \z\ 1+r }, 
Qz(fi{z),f2(z),e) = {{w,z) : \(w,z)\<e, fi{z)<w<f 2 {z)}, 
Q w (fi(w), f 2 (w),e) = {(w,z) : \(w, z)\ < e, fi(w) < z < f 2 (w)}, 

where we assume h{z) = 0(z 1+r ), f x (w) = 0(w 1+r ), h{z) - f 2 (z) = 0(z 1+ri ) } f\{w) - 
f 2 (w) = 0(w 1+ri ) for some r > and r x > r. In line with the definitions, the areas of the 
sets Q z (fi(z), f 2 ( z ), e) and Q w (fi(w), f 2 (w),e) are 0(e ri+2 ). Examples of the sets Q(r,e), 
Qz(fi(z), f 2 (z),e) and Q w (fi(w), f 2 (w), e) are shown in FigureEl By Theorem 14.11 below, if 
the stability index satisfies aj > —oo, then either the complement to £>f in B e (0) C Hj is 
empty, or it is the union of the sets Q z (fi(z), f 2 (z), e), or the union of Q w (fi(w), f 2 (w),e). 
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Figure 2: Examples of sets (a) Q(r, e), (b) Q z (fi(z), hiz), e) and (c) and Q w (fi(w), f2(w), e) 
used in the discussion of stability of Type A cycles. The sets are shaded. 



Whether it is the union of the sets Q z or Q w , depends on the difference between the con- 
tracting and the transverse eigenvalues. 

Examples of the sets B 9 S are shown in Figure [3] for different signs of bj and a,j — bj — 1. It 
can be observed that the sets are invariant with respect to the symmetry (w, z) — > (—w, —z). 
This is so, because the linearised systems f lTUj) evidently possess the symmetry, and the global 
maps (f)j are symmetric (being linear). 

For calculation of stability indices, we introduce the collection of functions {hij(y)} for 
1 < j < m, I < j, which are defined as followfl 

if ai — bi < 
if < ai - bi < 1 
if a\ — bi > 1 

The next theorem is the main result for Type A cycles, namely it gives the stability indices 
(jj for the collection of maps gj related to Type A cycles. Recall that the coefficients cij and 
bj of the map gj are related to the eigenvalues of linearisation of (CQ) near £j as ctj — / Gj 
and bj = —tjjtj. Recall, that Cj > and tj > for all j, therefore aj > 0. Following [16j, 
we denote 

Pj = min(aj, 1 + bj), 
P = Pi • • ■ Pm, and note that generically the non- degeneracy conditions 

aj^l + bj, bj^-1, p^l (19) 

4 If an index takes values 1, . . . , m, then the index value modulo m is understood here and below. Note 
that I in hij can take negative values. 



n jjKV) 



aihi + i tj (y) -ai + 1 

ai - k 
aihi +x Ay) - k 
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Figure 3: Examples of sets B 9 S (shaded) for Type A cycles for the cases (a) b± < 0, a\ — b\ < 1 
and (b) b x > 0, a\ — 6 X > 1. 

are satisfied. The Theorem below is stated and proved more precisely as Theorem IA.ll 

Theorem 4.1 For the collection of maps gj associated with a Type A cycle, the stability 
indices are: 

(a) If p > 1 and bj > for all j then (Xj >+ = oo and Gj- = for any j. 

(b) If p > 1, bj > — 1 for all j and bj < for j = Ji, . . . , Jl then Oj- = and 0"^+ are: 



a 



,-. = min h~- ( — — j — 1, where j = I 
s=j u ...,j l J ' s \ b s J { 



j, j < s 
j -m, j > s 



(c) If p < 1 or there exists j such that bj < —1 then o~j i+ = 0, = oo and the cycle is 
not an attractor. 

Proof: (a) Since p > 1, there exists a q > such that 

P = Y[(p 3 -q)>l. (20) 

By Lemma \A.2\ for sufficiently small e 

\gj(w,z)\ < \(w,z)\ p *- q for any (w,z), \(w,z)\ < e. (21) 

Consider (w, z) G H^. Inequality ( 12 ip implies that for a given 5, we can find an e > such 
that 

\9j,k(w,z)\ = \g jtl (w,z)\\g k (w,z)\ < \g jA (w, z)\\(w, z)\ p < 5 
for all 1 < j < m, k > 0, and \(w, z)\ < e. 



Therefore o"i + = oo and a\- = 0. The proof for j > 1 is similar. 

(b) Consider s = Ji for some /, whereby — 1 < b s < 0. For any small q > and 5 > we 
can find e, such that 

(w, z) G Q s ' s (e) = Q t (-\ z \-W»+i, \z\- l ' b °+\ e) 

implies 

\g s (w,z)\ > 5 

and 

(w, z) G H<£^ \ Q s ' s (e), where Q s ' s (e) = Q z (-\z\- l ^- q , \z\~ l ^- q , e) 

implies 

\g s (w,z)\ < 5. 

For simplicity, we ignore small q and say that if 

(w, z) G Q s ' s (e) = Q z (-\z\~ 1/bs , |^|~ 1/6s , e), 

then 

\g s (w,z)\ > 5, 

and if 

(w,z)EH^\Q s > s (e), 

then 

\g s (w,z)\ < 5. 

We also assume that e is sufficiently small and all estimates for (w, z) required in all the 
applied lemmas are satisfied. 

Denote by Q s ~ M (e) c the preimage of Q s ' s (e) under the map g s -i, and by Q^ s {t) C 
Hj 171 ^ the preimage of Q s,s (e) under g s o g s _ 1 o . . . o g^. By construction of the sets Q j,s (e), 
1 < j < s , for any (w, z) G Q^ s {e) the inequality 

|p* o ° • • - 9j(w,z)\ > 5 

is valid. 

The measure (area) of the set Q s ' s (e) can be estimated as 

£(Q s ' s (e)) = 0(e- 1/fes+1 ) = o(e^ s( - 1/6s)+1 ). 

By virtue of the definition of functions and due to Lemmas IA.3I and IA.41 the measure of 
the set <0 1,s (e) is 

£{Q l ' s {t)) = 0( e hiA-Ub s )+iy 
(Here and below, the measure of an empty set Q- ? ' s (e) is supposed to be e°°.) 
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Denote by Q 1_m,s (e) the preimage of the set Q 1,s ( e ) under a complete iteration g along the 
cycle, and by Q l ~ krri)S {e) the preimage under k iterations. The measure of the set Q 1_fcm ' s (e) 
is 

£/£l-fcm,«/ e N\ _ Q/ e /ii_ fcnilS (-l/6 e )+l\ 

Since p > 1, by the same arguments as used in the proof of part (a), if (w, z) 6 PI i? e 
does not belong to any Q 1 " fcm ' s ( e ) for all s — J±, . . . , Jl and k > 0, then 

|5f J: fe(w, z)\ < 5 for all 1 < j < m and > 0. 

By construction of the sets Q 1 ~ fcm ' s (e), if (w, z) e <5 1_A:m ' s (e), then 

(^■^(w, z)| > 5 for some j and k. 

By properties of the functions hij, the measure of the set Q 1,s (e) is larger than that of any 
other set Q 1 ~ km ' s (e) for k > 0. Since 

^IWjQ 1 ^)) = ( e min ^ 1 -(- 1 ^)+i) ) 

by definition of the stability index, for j = 1 the statement of the theorem, part (b), holds 
true. For other j the proof is similar. 

(c) Below we assume that 5 and e are sufficiently small, so that for < max(<5, e) 

the conditions of all lemmas to be applied hold true. 

We consider three following cases, which cover exhaustively all possibilities: 

• Suppose p < 1 and bj > —1 for all j. Since p < 1, there exists a q > such that 

m 

P = X\{p 3 +q)<l. (22) 

By Corollary IA.l( b). there exists r such that if (w,z) e Q(r, 5) then 

kfcKz)! > 5 (23) 

for some / and k. Consequently, if (w,z) G (g K )~ 1 Q(r, 6) for some K > 0, then the 
inequality (123]) is satisfied for k = k + K. The complement to Q(r, S) in B$ is the union 
of the sets Q z = Q z (\z\ 1+r , -\z\ 1+r , 5) and Q w = Q w (\w\ 1+r , ~M 1+r , <*)• Corollaries [O] 
and IA.3I imply existence of the limit sets 

Q l r = hm ^ and Qjf = lim {g k )~ l Q w ± 0, 

and by Lemmas IA.3I and IA.4] ff23l is satisfied for (w, z) e Ql! m U for some I and 
k. Therefore, <j\- = oo. 
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Figure 4: Examples of sets B 9 (shaded) for Type B or C cycles. 

• Suppose that b s < —1 for some s, and also at least one of the inequalities, p > 1, or 
a t — b t < for some t, is satisfied. Let the set Q s ' s (e) be defined as in the proof of part 
(b). Denote by R s ' s (e) the complement to the set Q s ' s (e) in B t , and by R 1 ~ km ' s (e) the 
preimage of R s,s (e) under g s o g s _ 1 o . . . o g\g k . By the arguments of the proof of part 
(b), the measure of the set R 1 ~ km ' s (e) is 

£(i? 1 ~ fcm,s (e)) = 0(e ,ll_fcm ' s( ' _6s ' )+1 ). 

By properties of the functions hij, lim^oo hi_ kmtS (—b s ) = oo, and therefore o\ - = oo. 

• Suppose p < 1, at — b t > for all t, and b s < —1 for some s. Let the sets R 1 ~ km ' s (e) 
be defined as in the previous paragraph. Since p < 1 and at — b t > for all t, the 
sets do not vanish in the limit k — > oo. By Corollaries IA.2I and \A.3\ the limit set 
i? hm = lmifc^oo R 1 ~ km ' s (e) does exist, and by Lemmas IA.3I and [A.41 inequality (|23|) is 
satisfied for all (w, z) e i? bm for some Z and /c. Therefore, a\- = oo. 

The proof for (Jj 7 -, j > 1, is similar. QED 



4.2 Type B and C cycles 

Examples of sets B 9 5 for Type B and C cycles are shown in Figure HI The sets are invariant 
with respect to the symmetries (w,z) —> (—w, z) and (w,z) — » (w, —z), because the linearised 
systems (fTUj) evidently possess these symmetries, and the global maps 4>j are symmetric due 
to linearity and invariance of the subspaces (w, 0) and (0,z). Therefore, we consider in this 
subsection only positive values of w and z; components of B 9 in other three quadrants are 
obtained on applying the symmetries. 
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As noted in [18], the maps gj : Hj — > Hj"j related to the cycles of Types B and C 
asymptotically to the lowest order have the form {Ew a \ Fw h] z) and (Ew bj z, Fw a: >), respec- 
tively. In the new coordinates {(,?]), ( = \nz and r] = Inw, the maps gj are linear: 

where the transition matrices of the maps are 

= \) I) 

for cycles of Types B and C, respectively. In the definition of stability indices, asymptotically 
small z and w (and therefore asymptotically large negative ( and 77) are assumed. Hence, 
we ignore finite \n.E and InF. 

Recall that the coefficients cij and bj of the matrices are related to the eigenvalues of 
linearisation of (p]) near as aj = Cj/ej and bj = —tjjtj. For the map g = g m o . . . o g 1 the 
transition matrix is M = M(g) = M m - • • M\. We introduce the notation: Mj t k and M^> 
denote transition matrices for the maps gj^ and g^\ respectively; M^ 1 ^ = .!//•••.! / ; : A{, 
M> v i = ( t 'ii)' i; i2) an d v 2 — (' u 2i5 t '22) denote eigenvalues and associated eigenvectors of the 
matrix respectively. If the eigenvalues are real, A{ > \ 3 2 is assumed. (Generically the 

eigenvalues are different.) 

A necessary condition for (w,z) to belong to B 9 S (see Subsection 13. 3j) is that g k (w,z) is 
bounded for all k. To leading order, the map g : (£, 77) — >■ ((, rj) is described by the transition 
matrix M(g). Due to linearity of the map, in the new coordinates the condition that the 
iterates (CkiVkY — M k ((,r]y are bounded by an S < (i.e., ( n < S and r) n < S, or, in the 
original coordinates, w < e s and z < e s ) generically is equivalent to lim^oo M k ((, rff = 
(—00, — oo)*. 

We denote 

U-°°(M) = {(x,y) : x < 0, y < 0, lim M n (x,yY = (-00, -oo)*}. 

n— >oo 

Evidently, C/ _00 (M) = implies B 9 S = 0. The conditions for U~°°(M) ^ are given in 
Lemma IB. II in terms of eigenvalues and eigenvectors of matrix M. They are: 

(i) the eigenvalues are real; 

(ii) Xi > 1; 

(iii) Ai > |A 2 |; 

(iv) v n vi2 > 0. 
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In terms of entries of a 2 x 2 matrix M = (a^) (where we assume an > 022) the conditions 
are (Lemma IB. 2D : 



0) 


(an - a 22 ) 2 

h ai 2 a 2 i > 


(24) 


(ii) 


fa n + a 22 . . \ , 
max 1 , a u + a 22 - a n a 22 + 012^21 1 > 1 


(25) 


(iii) 


an + 022 ^ „ 
2 > 


(26) 


(iv) 


a 2 i > 0. 


(27) 



For calculation of stability indices we introduce the sets Ur = {((, 77) \ ma.x((,r]) < R} 
and 

U R (a 1 ,(3 1 ,q 1 ;a 2 ,f3 2 ,q 2 ) = 
{{(, rj) e U R : («i + gi)C + {Pi + qi)v < 0, (a 2 + g 2 )C + (^2 + g 2 )^ < 0}, 

where R < 0. If all entries of the matrices Mj, 1 < j 1 < m, are non-negative, the stability 
indices of the related map can be calculated using the following Theorem which is stated 
and proved as Theorem IB. II 

Theorem 4.2 Let g be a map related to simple heteroclinic cycle of Types B or C and Mj, 
1 < j < m, its transition matrices. Suppose that for all j , 1 < j < m, all entries of the 
matrices are non-negative. Then: 

(a) If the transition matrix M = M m ■ ■ ■ Mi satisfies condition (ii) (see l[25\) ). then cr J>+ = 
00 and aj_ = for all j and moreover the cycle is asymptotically stable. 

(b) Otherwise, o~j 7+ = and cr,- _ = 00 for all j and the cycle is not an attractor. 

For calculation of stability indices of matrices with negative entries we introduce the 
following functions 



00, a > 0, > 0, 
0, a < 0, p < 0, 



P 



P 



1, a < 0, P > 0, — < -1, 

a a 

0, a < 0, P > 0, - > -1, 
a 

a a 

-- - 1, a > 0, p < 0, - < -1, 

0, a > 0, P < 0, ^ > -1 
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/ {a,P) = / + (-«, ~P) and 

r^{a,P) = f + {a,P)-f-{a,P) 

and prove in the following theorem that the set Bg(g) in H x (and similarly for Hj with 
j > 1) in the ((, rj) coordinates is U R (ai, P±, 0; a 2 , P 2 , 0). We denote the latter set in original 
coordinates (w, z) as 

U R (a 1 ,p 1 ,0;a 2 ,p 2 ,0) 

and then by Definition [5] 

r ln(£(U R \ U R {a x , ft, 0; a 2 , P 2 , 0))) - \n(£(U R (a u p 1} 0; a 2 , P 2 , 0))) 

a = hm = 

fl^-oc ln(£(^)) 

= min(/ index (a 1 ,/3 1 ),/ indcx (a 2 ,/3 2 )). 
The Theorem below is stated and proved as Theorem IB. 21 

Theorem 4.3 Let X be a simple heteroclinic cycle of Type B or C and Mj, 1 < j < m the 
associated transition matrices. We denote by j = j±, . . .ji the indices, for which some of the 
entries of Mj are negative; they are all non-negative for all remaining j . 

(a) If for at least one of j — + 1 the matrix does not satisfy conditions (i)-(iv) of 
Lemma \B.Bf then the cycle is repelling and Oj = — oo for all j. 

(b) If the matrices satisfy conditions (i)-(iv) of Lemma \B.2\ for all j = ji + 1, then 
there exist numbers Pi, a 2 , P 2 ), 1 < j < m, such that 

(i) U (a{, P(, 0; a\, P$, 0) ^ 0, 1 < j < m. 

(ii) For any S < and q > there exists R < such that 

M^\M^) k {U R {oP x , P{, -q; ai Pi -q)) C U s for all I, \<l < m, k>0, 

(Hi) 

hm(M^))(M^) fe (C,r/) t = (-oo, -oo) for all (Cv) e U (a{, P{, 0; o? 2 , p{, 0). 

fc— >oo 

(iv) 

U (a{,P{,0;ai,Pi,Q) = U- OD (M^)n f| (M^ l '^)~ l U fl f| (M^+ m ^)- 1 f/ . 

1<1<L 1<1<L 

(v) lf\ 2 >0 then 

U (a{, Pi, 0; ai, Pi, 0) = U~°°(M®) H f| {M^y x U Q . 

KKL 
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and the cycle is a Milnor attractor. 



Note that 



where a^' s ' are entries of the matrix M^' s \ 



Theorems 14. 1H4. 31 imply the following Corollary: 



Corollary 4.1 For simple heteroclinic cycles in M 4 
CTj = — oo for all j . 



(jj = — oo for some j if and only if 



4.2.1 Calculation of stability indices for Type B cycles 

Types B± and B± The cycles of Types and B^ have transition matrices 




Corollary IB.ll implies that if a < 1 or b < 0, then the cycles are not attractors and the 
stability index is — oo, otherwise they are attracting and the stability index is oo. 

Types B2 For cycles of Type B^ the product of transition matrices is 



with eigenvalues aia 2 and 1, and the associated eigenvectors (aia 2 — 1, &ia 2 + b 2 ) and (0, 1), 
respectively (for M 2 M 1; simply swap the indices 1 and 2 in the expressions to obtain the 
corresponding eigenvectors) . 

Theorems 14.21 and 14.31 imply to obtain the following classification: 

• If bi < and b 2 < 0, then the cycle is not an attractor and all stability indices are 
—00. 

• Suppose b\ > and b 2 > 0. 

— If a\a 2 < 1, then the cycle is not an attractor and the stability indices are —00. 

— If a,\a 2 > 1, then the cycle is locally attracting and the stability indices are 00. 

• Suppose b\ < and b 2 > 0. 

— If a\a 2 < 1 or b\a 2 + b 2 < 0, then the cycle is not an attractor and the stability 
indices are —00. 

— If a\a 2 > 1 and b\a 2 + b 2 > 0, then the stability indices are a% = / mdex (&i, 1) and 
a 2 = 00. 
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Type B 3 For cycles of Type B 3 the product of transition matrices is 

M 3 M 2 M 1= ( ^f 3 ). 

\ b 3 aia 2 + b 2 a x + b x I J 

Its eigenvalues are 010203 and 1 with associated eigenvectors 

(aia 2 a 3 - 1, b 3 a x a 2 + b 2 a x + 61), (0, 1). 

(For M 2 M X M 3 and MiM 3 M 2 , the quantities are obtained by cyclic permutation of the in- 
dices.) Theorems 14.21 and 14.31 imply obtain the following classification: 

• If b\ < 0, b 2 < and 63 < 0, then the cycle is not an attractor and the stability indices 
are all —00. 

• Suppose b\ > 0, b 2 > and 63 > 0. 

— If a±a 2 a 3 < 1, then the cycle is not an attractor and the stability indices are —00. 

— If a\a 2 a 3 > 1, then the cycle is locally attracting and the stability indices are 00. 

• Suppose bi < 0, b 2 > and 63 > 0. 

— If 010203 < 1 or 610203 + 6302 + b 2 < 0, then the cycle not an attractor and the 
stability indices are —00. 

— If aia 2 03 > 1 and 610203 + 6302 + 62 > 0, then the stability indices are 0\ = 
/ index (6i, 1), a 2 = 00 and a 3 = f ndcx (6 3 + M3, 1). 

• Suppose 61 < 0, 6 2 < and 63 > 0. 

— If aia 2 a 3 < 1 or b 2 a\a 3 + 6ia 3 + 63 < or 6ia 2 o 3 + a 2 6 3 + 6 2 < 0, then the cycle is 
not an attractor and the stability indices are —00. 

— If 010203 > 1, 620103 + 6103 + 63 > and 610203 + 0263 + 62 > 0, then the 
stability indices are a x = min(f ndex (6i, 1), / indcx (6i + 6 2 oi, 1)), o 2 = f ndcx (6 2 , 1) 
and ( T 3 = / index (6 3 + 6io 3 ,l). 

4.2.2 Calculation of stability indices for Type C cycles 
Type Ci Cycles of Type C-f have a transition matrix of the form 

"-(-)■ 

Corollary IB . 1 1 implies that the cycle is attracting and the stability index is 00 whenever 6 > 
and a + 6 > 1; otherwise it is not an attractor and the stability index is —00. 
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Type C 2 The product of transition matrices for cycles of Type C 2 is 

b]b 2 + a 2 h 



M t M 2 - , 

aib 2 a\ 

Denote by Ai and A2 eigenvalues of the matrix M\M 2 (which will be the same as those for 
MjM.). * 

Theorems 14.21 and 14.31 imply the following: 

• If bi < and b 2 < 0, then the cycle is not an attractor and the stability indices are 
—00. 

• Suppose that b x > and b 2 > 0. 

- If 

max(6i&2 + a 2 + a±, 2(6162 + a 2 + ai — a\a 2 )) < 2 
then the cycle is not an attractor and the stability indices are — 00. 

— Otherwise the cycle is locally attracting and the stability indices are 00. 

• Suppose bi < and b 2 > 0. 

- If 

(6i6 2 + a 2 — a{f + 4ai6i6 2 < 0, 

or 

max(&i&2 +02 + 0!, 2(6162 + a 2 + ai - oia 2 )) < 2, 

or 

6162 - ai + a 2 < 0, 
then the cycle is not an attractor and the stability indices are —00. 

— If none of the listed conditions are satisfied, the stability indices are 

a, = f ndcx ( M 2 + ax-A 2 aM ^ = /index( A 2 -M 2 -o 2 ^ 

62 61 

Type The transition matrix for cycles of Type is 

M = M 4 ' 1 = M4M3M2M1 = 
(6i6 2 + ai)(6 3 6 4 + a 3 ) + 6ia 2 6 4 (6 3 6 4 + a 3 )6 2 + 6 4 o 2 
a 4 6 3 (6i6 2 + ax) + a 2 a 4 6i 04^2^3 + a 2«4 

Denote by A 4 and A 2 eigenvalues and by vf' 1 and V2' 1 the associated eigenvectors of the 
matrix. The trace and determinant are: 

tr(M) = 6x62636/! + bib 2 a^ + 6 3 6 4 Oi + 6i6 4 a 2 + a 4 6 2 6 3 + aia 3 + a 2 a±, det(M) = a 1 a2a 3 a 4 . 

Theorems 14.21 and 14.31 imply the following: 
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If bj < for all j, then the cycle is not an attractor and the stability indices are — oo. 
Suppose that bj > for all j. 

- If 

max (trM, 2 trM - 2 det M) < 2, (29) 
then the cycle is not an attractor and the stability indices are — oo. 

— Otherwise the cycle is locally attracting and the stability indices are oo. 

Suppose that b\ < and bj > for 2 < j < 4. 

- If (EHD holds, or 

(trM) 2 -4detM < (30) 

or 

vgvg < (31) 
then the cycle is not an attractor and the stability indices are — oo. 

— If none of the listed conditions are satisfied, the stability indices are 

Oj = f indcx {v£ 3J /h j+3 ' j , -v£ 3J /h j+3 ' j ) where 

h j+3,j = ^+ 3 >V 22 +3J ' - v{+ 3 ' j vi+ 3 ' j and v{+ 3 ' j > 0, v{+ 3 ' j > 0. 

(v(i 3 ' : ' > and v\^ 3 ^ > can be assumed, because v x { v 12 2 > implies that 
v{l 3 ' j v{t 3 ' j > for all j). 

Suppose that b\ < 0, b 2 < 0, b 3 > and 6 4 > 0. 

- If (E9D, or flUD, or (J3IJI, or 

v 2 A 3 v 2 { 3 < 0, (32) 
holds then the cycle is not an attractor and the stability indices are — oo. 

— If none of the listed conditions are satisfied then the stability indices are 

a x = min(r dex (4 1 /^ 1 , -v^/fr 1 ), f mdcx (hb 2 + a u h), / index (6 l5 1)), 

a 2 = minir^iv^/h 1 ' 2 , -vg/h*), / index (& 2 , 1)), 
a 3 = min(/ indc > 2 2 2 7/i 2 ' 3 , -v 2 2 f/h 2 ' 3 ), / index (6i6 3 &4 + ha 3 + b 3 a,, hb, + a 4 )), 
a 4 = min(/ indcx (4 4 // i 3 ' 4 , -v%/h*>% 
r dex (hh + a 4) 6i). r de *(bib 2 b4 + b 2 a A + ai 6 4 , 6i& 2 + a a )). 
Suppose that b\ < 0, b 2 > 0, 63 < and 6 4 > 0. 



2s 




(33) 



holds then the cycle is not an attractor and the stability indices are — oo. 
If none of the listed conditions are satisfied then the stability indices are 




i+3,j 



J+3J 



/h j+3 ' j ). 



• Suppose that b\ < 0, b 2 < 0, b 3 < and 64 > 0. 

— If at least one of (I29p - (13"3"]) is satisfied, then the cycle is not an attractor and the 
stability indices are —00. 

— If none of the listed conditions are satisfied, the stability indices are 



a, = min(/ index (^ 1 //i 4 ' 1 , -vg/h** 1 ), r dex (6 1; 1), / index (6 x 6 2 + a x , b 2 ), 

r^ihhh + ai b 3 + ha 2 , b 2 b 3 + a 2 )), 
a 2 = min(/ index (4 2 // i 1 - 2 , -vtf/h 1 ' 2 ), f ndex (M 3 + a 2 , b 3 ), / index (6 2 , 1)), 
a 3 = min(/ indcx (4 3 /^ 2 ' 3 , -$?/h 2 ' 3 ), r^ihhh+b^+ha,, hh+a 4 ), / indcx (6 3 , 1)), 
a 4 = min(/ indcx (4 4 /^ 3 ' 4 , -^i/h z '% / index (fei6 4 +a 4 , h), / index (6 1 6 4 a 2 +a 2 a 4 , M 2 )). 



4.3 Comparison with earlier results 

Asymptotic stability of heteroclinic cycles has been previously examined in a number of 
papers. Type A heteroclinic cycles were considered by Krupa and Melbourne [TBI E] • I n 
the first paper [16J, cycles with negative transverse eigenvalues were investigated and the 
condition p > 1 was found to be necessary and sufficient for asymptotic stability of cycles. 
In the second paper [17] it was shown that cycles with some positive transverse eigenvalues 
are essentially asymptotically stable, if and only if p > 1 and the condition tj > —1 is 
satisfied for all j. This result is a special cases of our Theorem 14.11 The stability of Types 
B and C cycles with positive transverse eigenvalues was studied in [IB]. The conditions for 
asymptotic stability presented in this paper are equivalent to ours as given in subsections 
14.2.11 and 14.2.21 for the special cases bj > 0. 

A cycle of Type C 2 with one positive and one negative transverse eigenvalues was 
considered in [21]. Conditions (12)-(14) in [21], under which the set of points satisfying 
linifc^oo g k (w, z) = (0,0) is non-empty, are equivalent to our conditions that 0\ > —00 for 
b\ < and b 2 > 0. Existence of the set U (a{, (3{, 0; a J 2 , f3^, 0) was also noted in |21j. It was 
termed and was defined by the condition z > w l3+ , which we express as (3 + ( — i] < 0. 
This inequality implies that the stability index <j\ is equal to / index (/3 + , — 1). Substituting 
in it /3 + defined in the beginning of [21], Theorem 3.2], we obtain the value a 2 given in ff28l) 
(subject to the appropriate change of indices of tj and Cj). 
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5 Discussion 



Although it is natural to investigate cusp-like basins of attraction for heteroclinic cycles, 
to our knowledge this is the first paper to identify the algebraic order of the cusp as being 
an invariant of the dynamics — we use the stability index a(x) to characterize the local 
geometry of basins of attraction near invariant sets X in general, and heteroclinic cycles in 
particular. This quantity might be especially useful in describing the local structure of a 
range of invariant sets, for example, for riddled and intermingled basins of attraction [fl [6]. 

In the latter part of the paper we have calculated how the stability indices depend on 
the cycle structure and eigenvalues for simple (robust heteroclinic) cycles in K 4 . Clearly, 
transition matrices can be used to study the stability of simple cycles in higher-dimensional 
systems or for more complex cycles; however, we expect such a classification to be so complex 
that the results can hardly be enlightening — and so we have not attempted this. 

Our approach should give some insight into the structure of heteroclinic networks [3j 18] 
and heteroclinic switching [5j [TOj [121 1131 EH]- For some cycles (of Type A where at least 
one difference between transverse and contracting eigenvalues is negative, or of Type B 
where at most one transverse eigenvalue is positive) we find one of the stability indices 
to be oo. Consequently, if a heteroclinic network involves such a cycle, no switching is 
possible for perturbations on that particular connecting trajectory (almost all trajectories 
that are near a connection where o = oo will stay near the cycle for all t > 0). Some of 
the results in this paper may be extended to general robust heteroclinic cycles; it seems 
plausible that Corollary 14. II holds for simple heteroclinic cycles in IR n where all connections 
are one-dimensional manifolds, because for such cycles the Poincare map along the cycle 
becomes linear after the same change of coordinates. It should also be possible to extend to 
cases of compact but not finite symmetry, though again this is likely to be quite involved 
owing to the complexity of heteroclinic cycles between relative equilibria. 

Finally, we emphasise that there is no a priori reason why the limits in the definition of 
the stability index exist. We give below an example where the stability index can be shown 
not to converge. This is probably not a generic situation though; the example below is 
highly degenerate, and the generic conditions for the cycles in R 4 , as detailed in the previous 
section, all result in computable stability indices. 

An example where <j + and cr_ do not exist. We consider a (non-invertible) map 
M : [0, oo) — > [0, oo) with X = {0} and consider its basin of attraction N = B(X). Define a 
sequence = exp(— 2 k ) and 

if e 2 fc +2 < \y\ < e2fc+i 
if e 2 fc+i < \y\ < e 2fe 
if y = 

Then 

v t2k — C2fc+1 _ 1 

^e 2k > — 1 — e 2 fc, 

^2k 



M{y) 



(e 2 k~i(y - e 2 fc+2) 



€2k(y — e2fc+l))/( e 2fc+l — £2k+2), 

o, 

0. 
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whereas 

^e 2k+1 < - £2fc+l- 

Hence, ln(£ £2 J/ ln(e 2fe ) < ln(l - e 2fc )/ ln(e 2 fe) and ln(E £2fc+1 )/ ln(e 2fc+ i) > 1, and therefore the 
limit defining cx_(0) does not exist; it can similarly be shown that cr + (0) does not exist. This 
example can clearly be extended to a continuous map M. Although it is not easy to see how 
to extend to a map that is differentiable at y = with the same properties, it may well be 
possible to produce a smooth example in dimension two or more. 



Other global measures of stability The stability index a(x) describes the local geome- 
try of the basin of attraction of an invariant set X. One can define global and local stability 
numbers of a flow invariant set X as follows: 

n (X)-Km emX)nB{X)) 
n glob (X)-hm £mx)) , 

l(B e (X)nB 5 (X)) 
^ioc(A) = hmhm - t— ■ , 

where Bs is defined as in ([8]). Clearly, from the definition one can verify that < n g i b(X) < 
1, < n\ oc (X) < 1, and n g i Q b(X) = 1 if and only if X is p.a.s. for a local basin of attraction. 
Note however that stability number is not an invariant of the dynamics; we believe that 
only the classification into whether n(X) = 0, < n(X) < 1 or n(X) = 1 and scaling 
properties will be invariant under smooth conjugation. For a heteroclinic cycle comprised of 
one-dimensional connections, the stability number can be related to its stability indices by 
the following: 

nglob(x) - efW) ' 

where C l {sj) is the length of the connection sj. 
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A Type A cycles 

In this Appendix we present a proof of the main theorem for calculation of stability indices 
for Type A cycles. Consider the map g = g m o . . . o g x : M? — > IR 2 , where^] 

gj(w,z) = {gj{w,z),g z j {w,z)) = (A jW a ^\z\ a ^ + Bj\w\ bl * z* 2 * , C jW a ^\z\ a ^ + Dj\w\ b ^ z^). 

(34) 

For generic cycles of Type A, AjBjCjDj ^ and AjDj ^ BjCj holds true for all j. In 
the first two lemmas in this appendix we assume that aij + a 2 j < bij + b 2 j. No generality 
is lost, because the expression f fTB"]) for the map gj(w,z) is invariant under the transforma- 
tion (Aj, Bj, Cj, Dj] aij, a 2 j, bij, b 2 j] w, z) — > (Bj, Aj, Dj, Cj; b 2 j, by, a 2 j, a^; z, w). Recall that 
gi,k = gi o . . . o gi o g k . 
Let us define 

Q{r,e) = {(w,z) : \{w,z)\ < e, \z\ > \w\ 1+r , \w\ > \z\ 1+r }, 
Q z (fi(z),f 2 (z),e) = {(w,z) :\(w,z)\<e, f x (z) <w < f 2 (z)}, 
Qw(fi(w), f 2 (w),e) = {(w,z) : \(w,z)\ < e, f^w) < z < fe(w)}, 

where we assume h(z) = 0(z 1+r ), f^w) = 0(u> 1+r ), f x {z) - f 2 (z) = 0(z 1+ri ), fi(w) - 
f 2 (w) = 0(w 1+ri ) for some r > and r x > r. 

We start the study of stability by proving a few lemmas about properties of the maps gj. 

5 By virtue of (|13|) . a,2j — and — 1. However, the first two lemmas in this Appendix are proved for 
arbitrary and bij . 
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Lemma A.l Suppose ay + a 2 j > 0. For any q > and r > satisfying 

(1 - r)(6y - ay) + 6 2 j - a 2 j > 0, 6y - ay + (1 - r)(6 2 j - a 2 j) > 0, ay + a 2j - - g > 
and 



r < mm 



Q 



2a±j 



Q 



2a 



there exists an e 7 - > ; such that 



■2.1 



aij+a,2j —g\ 



and 



9j(Q(r,tj)) C Q{r,e j 
\gj{w,z)\ < \{w,z)\ ai ^- q for any (w,z) G Q{r, e,) 

\gj(w,z)\ > \(w,z)\ ai3+a2i+q for any (w,z) e Q{r,€j) 



(35) 

(36) 
(37) 

(38) 



Proof: At least one of the coefficients, ay or a 2 j, is positive. Assume ay > 0. The inequality 
ay + &2j < bij — b 2 j implies 6y — ay + b 2 j — a 2 j > and hence at least one of two differences, 
6y — ay or b 2 j — a 2 j, is positive. For the sake of definiteness we assume without loss of 
generality that 6y — ay > 0. Set tj < 1 satisfying the following inequalities: 



€ j 



(1— r)(hj— aij)+&2j— a>2j 



< mm 



2B j 



2D, 



ef 1(3^/2, 3^/2)1 <1, 



r(ai j +a 2 j-g) 



efK< 1(^/2,^/2)1, 

(to,z) G Q(r,ej-) 
"' + /,' ; V". "''.r', + DjW hli ~ aii z h%i ~ a2j ] 



A; 



3d 



and 



whera^if = 2 a+b+q+1 . 
Assume that 

and re-write (TT81) as 

gj(w,z) = w aij \z 
Due to (|32J, 



(39) 

(40) 
(41) 
(42) 

(43) 



\w 



bi 3 -aij z b 2j -a 2 j I < I ^ I (&ij -ay )/(l+r)+o 2j -a<y < e ( 1 - r )( 6 li- a lj)+ 6 2j-a2j 



< mm 



-4; 



2R- 



a, 



2D, 



(44) 



Here, and below, we use the norm \(w, z)\ = (w 2 + z 2 ) 1/2 - If a different norm is employed, the proofs 
remain similar but some constants will be modified. 
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Therefore, due to (|35]h (|40|) and (|43 



\gj{w,z)\ < \w ai 'z a2 ^3A J /2,3C J /2)\ < |2| a ^' /(1+r)+a ^ 1(3^/2, 3^/2)1 
< \(w,z)\ ai ^ 1 - r)+a ^- q/2 e q / 2 \(3A j /2,3C J /2)\ < \{w, z )\ ai ^- q . 



(45) 



Thus, ( 137|) is proved. 

The inequalities ( 144"|) . ( )4"Tj) and ( l4"5j) imply that for any g > and r > the €j can be 
chosen so that 



gf(w,z) 



> 



A; 



3Cj 



> e 



■ {a ^- g) >\9 z Aw,z)\ r . 



Hence 

and similarly 



\ 9 y(w,z)\>\g*(w,z)\ 



l+r 



\g*(w,z)\>\gj(w,z)\ 1+r . 

Therefore f )36|) holds because of (T45|). 

Combining (141 |) . P?]) . ( 155) and (132 ]) we obtain 



|s>,z)| > 1^^(^/2,^/2)1 > |^r«( 1+ ^|(A,/2,^/2)| 
> |^| a «( 1 + r ')+^+3/2 e -'?/2|(^y2 )C ' j /2)| > urizl'w+^+a. 

Assume that a 2j < 0. Estimates (111 ) . ( B3"]) . p)]) and (H2 ]l imply that 

\ gj (w,z)\ > \w^z^(A j /2,C j /2)\ > Iwr^/^UjfrCjWl 
> | ^ ^1,+^ (i-r-) | | > |u;| a ^ +a ^ (1 - r)+9/2 eT' ? / 2 |(^/2,C J /2)| 



> K\w 



aij+a,2j+q 



If a2j > then 



\gj{w,z)\ > \w\ ai i +a2 ^ 1+r) \(A j /2,C j /2)\ > K\w\ ai ^ +q . 



Thus, 



\gj(w,z)\ > -K(\z\ a v +a v +q + \ w \ ai i +a v +q ) > \(w,z)\ ai ^ +q . 



(46) 



(47) 



QED 



Corollary A.l Suppose that the conditions of Lemma 14.11 are satisfied for r > 0, q > 
and for all 1 < j < m 

( a ) V ni<j< m ( ai i + °2j ~~ l) > ^ ^ en f or anv ^ > there exists an e > s«c/i that 
\gi,k{w, z)\ < 5 for all (w, z) G Q(r, e) and for any 1 < I < m and k > 0. 
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(b) Denote 8 = mini<j< m ej. If rii<j<m( a ii + °2i + O) < 1 then for any 5 < 5o and any 
(w,z) G Q(r, 5) there exist I and k such that \gi^(w, z)\ > 5. 

Lemma A. 2 Suppose aij, a 2 j, bij, b 2 j > 0. For any q > there exists an e > such that 

\gj(w,z)\ < \(w,z)\ aiJ+a2j ~ q for any (w,z) with \(w,z)\ < e. (48) 

Proof: Let e < min((|A j | + \Cj\ + \Bj\ + {Djl)" 1 ^, 1). Thus, for such (w,z) we have 

\9j(w,z)\ < (\Aj\ + \Cj\)\w a ^^\ + (\Bj\ + \Dj\)\w b ^z b v\ 

< (\Aj\ + \Cj\ + \Bj\ + \Dj\)\(w,z)\ ai i +a2 i < \{w,z)\ a ^ +a ^- q . 

QED 

The following lemmas and the main theorem of this subsection consider the special case 
relevant to Type A cycles, where the maps gj(w,z) have a 2 j = and b 2 j = 1, i.e. they 
simplify to 

gj(w, z) = (A jW aj + Bj\w\ bi z, CjW a J + Dj\w\ bj z) (49) 
(we set d\j = CLj, b\j = bj and do not assume necessarily that aj < bj + 1). 

Lemma A. 3 Suppose aj — bj > I. 

(a) For any q> 0, fi(z), fziz), V andp, such that p > p > 1 and 

f[(z) = 0(^- 1 ), h{z) - f 2 (z) = 0(z% as z^ (50) 
there exist fi(w), f 2 (w) and ej > such that 

f[(w) = 0( w a *- b i- 1 ), AH - f 2 (w) = 0{w a ^) forw^ 0, 

and 

ft(^(/iH,/2H,e 3 )) C Q z (fi(z)J 2 (z),e^- q ), 
9j (QMiHJ2(w),e J )) D Q z (f l (z),f 2 (z),e a / +q ), 
\ 9j (w,z)\ < \(w,z)\ a ^, \ 9j (w,z)\ > \(w,z)\ a ^ 
for all (w,z) E Q w {h{w),f-2.{w),ej). 

(b) For any q > 0, f\{w), f'2(w), p and p, such that p > p > 1 and 

f[{w) = OK" 1 ), AH - f 2 (w) = OK), as w 0, (51) 
there exist fi(w), f 2 (w) and €j > such that 

f[(w) = O^ 4 '- 1 ), /iH - f 2 {w) = 0(w a ^) asw^O, 

and 

9j(.Qw(.fi(w),f 2 (w),ej)) C Q w (fi(w), f 2 (w),e a /~ q ), 
9j(Qw(fi(w),f 2 (w),ej)) D Q w (f 1 (w)J 2 (w),e a / +g ), 
\ 9j (w,z)\<\(w,z)\ a ^, \ 9j (w,z)\>\(w,z)\ a ^ 
for all (w,z) G Q w (fi(w), f 2 (w), e,-). 
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Proof: (a) Let f 1,2(1*}) be respectively the solutions of 



and 



AjW aj + B 3 \w\ b if 2 (w) = f2(C jW a i + D j \w\ b *f 2 (w)). 



(52) 
(53) 

(54) 



Substitution of w = w + 5w into (1521) implies 

Aj(w+5wp +B j \(w+5w)\ b ^ hiw+dw) = f^C^w+Sw)^ +D j \(w+6w)\ b * friw+Sw)). (55) 

Subtracting now (|52|) from (I55I) . dividing the result by 5w and taking the limit 5w — > 0, we 
obtain that 



The functions /i j2 (io) ar e defined for any small w and due to (150]) 



5; 



Subtraction of (1531) from (1521) yields 



A?'( a j &i)„„a,— 6,-1 _|_ q/ ajp-ftj—^ 



-w 



^M^M - / 2 H) = /i(c> a ' + £>p/i(«0) - /i(^ aj + £>n/ 2 H)+ 

A(c> aj + d>| 6 '/ 2 m) - / 2 (c> aj + d>| 6 '/ 2 M). 

Because of (J50~j) and ([5 



(56) 



fx{CjW ai + D>|*/iM) - fi(C jW ^ + Di\wpf 2 {w)) 



D,A„ 



(57) 



and 



A(C> ' + D>| 6 i/ 2 H) - / 2 (C>^ + £>>| b '/ 2 (it;)) = 0(^"). 
Hence (156|) takes the form 

B^ihiw) - h(w)) = oK^-^oaiH - AH) + o(^). 

Since p > 1 and aj > 1, the first term in the r.h.s. of this expression is asymptotically 
smaller than the l.h.s. and it can be ignored. Thus, 

/iH-/ 2 H = o(^). 

The asymptotic relation ( 1541) implies that there exist e and K such that 



B 3 



< K\w 



ajp-bj I 
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for (w,z) G Q w (fi(w), f 2 (w),e). Suppose that 6j satisfies 



e 3 < min(e, 1), e) (\CA + Mi + (\ Bj \ + \Dj\)k) < 1 



K(tj {p X) < min 





A, 




A. 


D j 






B j 



B, 



4 LI 
, e. 



at—bi—X 



< 



Ba 



A j 



Then 



\gj(w, z) | = \w a J (Aj + Bjw^-^z, Cj + Djw^'^z 



< M^IQI + 1^1 + (1^1 + |£>,-|)/0 < \(w,z)\ a >- q 



and 



\gj(w,z)\ > \w aj (Cj + DjW b >~ aj z)\ > \w\ a i 



c . _ DA 



' B 3 



> 2\w\ a > +q > \(w,z)\ a i 



dj+q 



From (1521) and (|53|) . the proof of part (a) is complete, 
(b) Set fi, 2 (w) to be respectively the solutions of 



and 



C jW a * + Djlwfihiw) = fi{A jW a i + B 3 \w\ b >h{w)) 



C 3 w a > + D 3 \w\ b if 2 {w) = h{A 3 w a i + BAw\ b *h(w)). 



The remainder of the proof is similar to 



QED 



Corollary A. 2 Lemma \A.3\ implies that the maps f] — > f] are monotonic in the following 
sense: 

(a) Let fo be such that fi(z) < fo(z) < f 2 (z) for all \z\ < e* 3 q . Define fo(w) by analogy 
with fi i2 (w) in part (a) of that lemma. 

# AM > f 2 (w), then h(w) > f (w) > f 2 (w), 

Iffi(w) < f 2 (w), then f^w) < f (w) < f 2 (w) 

for all \w\ < 6j. 

(b) Let /o be such that fi(w) < fo(w) < f 2 (w) for all \w\ < <t 3 q . Define fo(w) by analogy 
with fi t2 (w) in part (b) of that lemma. 

Iffi(w) > f 2 (w), then A(u>) > fo(w) > f 2 (w), 

If fi{w) < f 2 {w), then f^w) < f (w) < f 2 (w) 

for all \w\ < t 3 . 



38 



Lemma A. 4 Suppose that < a 3 : — b 3 < 1. 

(a) For any q > 0, f\{z), f 2 {z), p and p, such that p > p > 1 and 

f[(z) = O^ 1 ), h{z) - f 2 (z) = O(^), as z^ 0, 
i/iere exzsi fi(z), f 2 (z) and e 3 > snc/i that 

f[( z ) = 0(z 1/(a ^ b ^ 1 ), hiz) - f 2 (z) = 0{w^-^ +1)/ ^- h ^) forz^ 0, 

and 

gj {Q z (Uz),f 2 {z),ej)) c Q,(/i(*), / 2 (*), e^-^" 9 ), 

gsiQzifiWtMzUj)) D Q.(/i(4/ 2 (^ef (a ^ )+<? ), 

|<7>,z)| < |(«;,«)|^-^-«, |<?>,*)| > \(w,z)\ a ^- b ^, 
for all (w,z) E Q z (fi(z),f 2 (z),ej). 

(b) For any q > 0, fi(w), f 2 (w), p and p, such that p > p > 1 and 

f[{w) = OK 1 ), /iH - /» = OK), as w 0, 
i/iere exzsi fi(z), f 2 (z) and e 3 > snc/i £/ia£ 

f[( z ) = 0(z 1/(a ^-^ ) - 1 ), - / 2 (z) = 0(z (o ^-^ +1)/(o ^-^ ) ) /or z -> 0, 

and 

9i(Q*(fi{*),M*Ui) c Q^/iH,/^),^"^), 
9j(Q,(h(*)J*(*Uj) d ^(/iW,/2H )e f (ory+? ), 

|<?>,z)| < |(«;,«)|^-^-«, |<7>,4 > IK*)!^"^, 
/or a// G Q 2 (/i(2;),/2(2;),e j ). 

Proof: (a) Let fi, 2 (z) be solutions of 

V? (*) + = MCjffiz) + ^|/i| 6 ^) 

and 

A^OaO +^1/2^(4* = h{Cj 2 j {z) + D j \f 2 \ b ^(z)z). 
The functions fi, 2 (z) are defined for any small z. Note that 

/ D \ V(<»J— M 

/ x 2 (2) = - z V( a i- h i) + o(£ (ajP ~ aj+1)/(aj ~ 6j) ). 
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We substitute z — > z + 5z into ( 160]) . subtract ( 160]) from the obtained equation, divide the 
result by 5z and take the limit 5z — > 0. Since 



(1 j b j : ^ V -^4 -j 



f?(z+6z)-f?(z)=f?(z) 11 + 

and similar estimate holds true for 6,-, we obtain that 

1 



gl/(aj-bj)-l _|_ Q/ 2 ,(o 3 -p-Oj+l)/(oj-6j)-l> 



Subtracting ( IdTT) from ( 1601 we obtain 



(63) 



(64) 



Since 

ft 3 - P = (A - f2)o(z^- l ^-»>)), 
z(\h\ b > - |A| 6 o = (A - A)0(^- 1V ^), 

+ D,\f^z) - f 2 {Cj? + ^-lAl 6 ^) = 0(^/(^-6,)) 

and 

MCjf? + ^lAlN - f2(cj? + ^-l/ 2 | 6 ^) = (A - f 2 )o(z^-^- b ^), 

(164"]) implies that 

fi{z) - h{z) = o(z^- a ^ +l)/ ^- b ^). 

Due to ([62]) there exist e and K such that 



< tf|z|feP->i>/(«i-«y) for (to,*) G Q^z), f 2 (z),e). 



Suppose that ej satisfies 

ej < min(e, 1), A 



2Bj 



A,- 



aj/(aj-bj) 



|C,| + 



J 4 



2A- 



Oj/(oj-6j) 



< mm - 
2 



+ (1^1 + 1^1)^) <1, 
A,- 



Then 



>1 and eV^-V- 1 < 



A, 



< 



2B< 



A, 



\ gj {w, z)\ = \w a *{Aj + BjW*i- a 'z, Gj + DjW^-tz 

Oj/(Oj-6j) 



laj/(aj-bj) 



B, 



(I^-i + i^dk 



< w,2 



aj/(a,j-bj)-q 



(65) 
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and 



\gj(w,z)\ > (!■■'■(() • I ),(>■'■■ "•:) 
2 C '~~£ 



> 2 M 



aj / (a,j -bj)+q 



> (w,z) 



ctj/(aj-bj)+q 



Hence part (a) is proved. The proof for the part (b) is similar. 



(66) 



QED 



Corollary A. 3 Lemma \A.4\ implies that the maps f] — > fi are monotonic in the following 
sense: 



(a) Let f be such that f\(z) < fo(z) < f 2 (z) for all \z\ < e 



Define f (z) by 



analogy with fi )2 (z) in part (a) of that lemma. 

Vfi(z) > h{z), then fi(z) > f (z) > f 2 (z), 
Iffi(z) < f 2 (z), then f\(z) < f (z) < f 2 {z). 

for all \z\ < 6j. 

(b) Let f be such that f\{w) < fo(w) < f 2 (uj) for all \w\ < ey 1 b ^ q . Define fo(z) by 
analogy with fi l2 (z) in part (b) of that lemma. 

Ifh{z) > f 2 (z), then f x (z) > f Q (z) > f 2 (z), 

Iffi(z) < f 2 (z), then f\(z) < f Q (z) < f 2 (z). 

for all \z\ < 6j. 

Lemma A. 5 Suppose aj > and aj — bj < 0. Then for any eo > and r > there exists 
an e > such that 



gj(w,z) e Q(r, e ) for any (w,z) with \(w,z)\ < 



Proof: Let e satisfy 



< min 



A* 



2B j 



C; 



2D, 



6^1(3^/2,3^/2)1 <e 



and 



e raj < min 



-l-r 



A j ( 3Cj 



-l-r 
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Therefore, if |(tu, z)| < e then 

\gj(w,z)\ = \w aj \\(Aj + /,' ; V " + DjW bj ~ aj z)\ < e a *\ (3Aj/2, 3Cj/2)\ < e , 

\gj(w,z)\^/\g*(w,z)\ < l^l ^r^ < 1 

and similarly 

\g*(w,z)\ 1+r /\gr(w,z)\<l. 



QED 



Lemma A. 6 For any q > there exists an tj > such that 

\gj(w,z)\ < \{w,z)\h where = min(a,/2, \bj\q/2, |1 + &?|/2), 
for all (w,z) G B e , \Qi, where 

ft J0 

I Q*(-kr 1/6 *~ f ,N~ 1/6jH Uj) ifbj<0. 

Proof: If 

(l^'l + IQI)^ 72 < V2, (ISyl + \Dj\)# m < 1/2, (1^1 + |^.|)6^ l/2 < 1/2 
then one can verify that 

\ gj (w,z)\ < (\Aj\ + \C s \)\w\^ + (\BA + \D J \)\w b >z\ < \(w,z)f>. 



QED 



Lemma A. 7 For any q > 0, r > and e > there exists an €j > such that 

gj(w, z) e Q(r, e) for all (w, z) e B ej \ (Qi UQ 2 U Q 3 ), 
where Q\ is defined in Lemma \A.6[ 

3 if a,j — bj < 

Q 2 = <( Q*(-(ff) ai 2 ai - M a2 , + M" 2 , Ci) i/0 < a, - b, < 1 

Q4-^j-w a3 - H a4 ,-^u; Q3 + Iwl** 4 , tfcij-bj > 1 

i/ dj - bj < 

= <( ^H^V 1 - M a2 ,-(§) ai 2 Ql + M° 2 ,e,) z/0 < aj -b 3 < 1 
Qu,(§-^ Q3 - \w\ a \ %-w° 3 + M a Vj) i/aj- - 6 3 - > 1 



and 



1 a,(r + 1) - 2a, + 2 a 7 (r + l) 

a j ~ bj 2(aj — bj) 2 
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Proof: We start with the proof of existence of e w such that in the case aj — bj > 1 the 
condition 

{w, z) E B ew \ {Q 1 U Q 2 U Q 3 ) 

implies that 

\gf(w,z)\>\g](w,z)\ 1+r . (67) 
If \AjW a * + Bj\wpz\ < \AjW a */2\ then 

\gj(w,z)\ > |BX i(r+1)/2 l and < (l^'l + 

For | w %(^+i)/2| < I^KIQI + fiAjDj/Bjl)-*- 1 the inequality (JSZD holds true. 
If \AjW a * + Bj\wpz\ > \AjW a '/2\ then 

\gj(w,z)\ > min(|A j /2|, |Sj/3|)max(|tw^|, \w b *z\) 

and 

\glj(w,z)\ < 2max(|C i |, \Dj\) max(\w aj \, \w bj z\). 

Hence if 

\(w,z)\ rf, t <min(|A / /2US i /3|)(2max(|<7 i |,|£) i |))-'-- 1 , 

where (3j is defined in Lemma [A. 6 1 then (1671) holds. 

The proof of existence of e z such that \gj(w,z)\ > \gj (w , z)\ 1+r is similar. Denote by 
6j the €j from Lemma TA.6I Then £j = min(e w , e 2 , ej), by the definition of the set Q(r,e), 
satisfies the the condition of the lemma. For the case < aj — bj < 1 the proof is similar 
and is not presented, for the case aj — bj < the statement of the Lemma follows from 
Lemma \KM QED 

Let the collection of functions {hij(y)} for 1 < j < m, I < j, be defined as follow^: 
hjAv) = v> 

{oo if a>i — bi < 
aihi+i j(y) - ai + 1 
if < a, - k < 1 
a; - k 
aihi +ljj (y) -k if a t - b x > 1 

This collection has the following properties: 
(a) 

hjiVo + Vi) = hi d (y ) + aijyi, 
where ay = Yli< s <j m ax(a s , a s / (a s — b s )) is defined for I < j. 
(b) If there exist J such that aj — bj < 0, then hij(y) = oo for I < J. 



7 If an index takes values 1, . . . , m, then the index value modulo m is understood here and below. 
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(c) If %- mJ > 1, then hi- m j{y) > h d {y) for any y > 1. 

(d) If cij-mj > 1, then lim^-oo hij(y) = oo for any y > 1. 

The next theorem gives the main result for Type A cycles, namely it gives the stability 
indices <jj for the collection of maps gj related to Type A cycles. The coefficients aj and bj 
of the map gj are related to the eigenvalues of linearisation of ([T]) near £j as a, = Cj/ej and 
6-,- = —tj/ej. Recall, that Cj > and ej > for all j and therefore aj > 0. Following [16], we 
denote 

Pj = min(aj, 1 + bj), 

P = Pi ■ ■ ■ Pm, and note that generically the non- degeneracy conditions (|T9|) apply. 



Theorem A.l (reproduces Theorem \4-l\ ) For the collection of maps gj associated with a 
Type A cycle, the stability indices are: 

(a) If p > 1 and bj > for all j then o~j^ + = oo and o~j- = for any j. 

(b) If p > 1, bj > — 1 for all j and bj < for j — Ji, ■ ■ ■ , Jl then aj - = and &j i+ are: 

a ji+ = min hj, s ( -— 

s=Ji,...,J L \ O s 

(c) If p < 1 or there exists j such that bj < —1 then &j j+ = 0, aj - = oo and the cycle is 
not an attractor. 

Proof: (a) Since p > 1, there exists a q > such that 

H( Pj -q)>\. (68) 

By Lemma IA.21 for any j there exist €j such that 

\gj(w,z)\ < \(w, z)\ pj ~ q for any (w,z) with < e,-. (69) 

For a given 5, choose an e > satisfying 

j 

e Pi,j < m m(5, €j, 1) and pij = J^[(Ps — q), for all 1 < j < m. (70) 

s=l 

Consider (w,-2) G /Z^ . If | (w, z) | < e, then fl69l) and flTUj) imply that |(7 J)0 (w, z)| < 5. Due 
to fl68l) . < e and hence ^^(w, z)| < 5 for all < j < m — 1, A; > 0. Therefore 

o"i + = oo and ax _ = 0. The proof for j > 1 is similar. 



(b) First, let us prove that 



< h ljS [ ~ ) - 1 + q u (71) 
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Figure 5: Examples of the sets (a) Q s ' s (e s ), (b) Q s 1,s (e s _i) for C s _ii} s _i < and < 
a s -i - b s -i < 1 and (c) Q s ~ 2,s (e s - 2 ), for A s _ 2 B s - 2 < and a s _ 2 - b s - 2 > 1. 

where q~\ is any small number and s — J\ for some I. Denote qi = q\/ai tS . Assume that q 
satisfies (168|) . Define the sets Qi' s (ej) by the following rule: 

Q s ' s (e s ) = Q z {-\z\- l/b ° +q \ \z\- l ' b ^\e s ). 

For a given Q J ' +1 ' s ( e j+i) the set Qi> s (ej) is 

{0 if Oj - 6 3 - < 

Q z (fi(z), f 2 (z), e 3 ), e j = ^M^-^) if o < — 6, < 1 (72) 
Qw(fi{w), f 2 (w), €j), 6j = effi~ q ' if - > 1, 

where f\ and / 2 are the functions defined in Lemmas IA.3I and IA.41 Denote by ej the 6j from 
Lemmas I A. 31 and [ A .41 and set e = min^ e^ ai ' J ; where ay = ]J Ks<J (max(a s , a s /(a s — b s )) — q). 
Examples of the sets Qi ,s (e s ) are shown in Figure 

For any 5 > we can find an e > such that 

\g s {w,z)\ > 5 for all (w,z) G ^H^" 176 ^ 91 , M~ 1/6s+<zi , e)- 
Hence, if e < min(eo, e 1//a ^ 1 ' s - ) ) then 

\g s o ...o gx(w, z)\ > 5 for all (w, z) e Q 1,s (e). 

Since 

and hi )S (l/b s + qi) — h\ )S {l/b s ) = q~\ the inequality (I7T1) is proved. 
Second, we prove that 

a l>+ > h ltJ (-l/bj) -l-q u (73) 
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where q\ is any small number, g x = q\/a l j, and J is the value of J\ where 

min (/ii j 

is achieved. Assume that q satisfies (1681) and r and g satisfy the conditions of Lemma IA.1I 
for all j. Set 

Q s ' s (e s ) = Q z (-\z\- 1/b °- q \ |^|- 1/6s - 91 , e s ) for s = 1, . . . , J z , 

Q^' 1 (ei) = Qi where Q, I = 2,3, and ei are defined in Lemma IA.7I for j = 1. The sets Q^ s 
and Qf 1 , I = 2, 3, are defined by f[T2"|) . Denote 

{oo if a x — b\ < 

(ai(r + l)/2 - ax + l)/(ai - &i) if < a x - &i < 1 
ai(r+l)/2-6i if ai - &i > 1. 

Since y > 1, by property (d) of the functions hij(y), there exists k > such that hi_ mkjl (y ) > 
h 1 j(—l/bj). Hence by Lemmas IA.6I and IA. 71 

if (w, z)eB e \ ((U < s < fc , 1 < l < L Q 1 - sm > J <) U Q\- km < 1 U Q^" 1 ' 1 ) then g k (w, z) G Q(r, e), 
where e kp < e, (3 = ]\i<j< m Pj- Since 

£(Q 1 - J (e)) = 0(e^ ( - 1/b ^ l)+1 ), 
^(Q^ fem,1 (e)) = 0(e hl — m(ko)+1) = o ^^ J -(-i/6 i - gi )+i) 

and 

£(g 1 " sm ^(e)) = 0(6^— .Jf(-V6jr-9i)+i) = ( e ^i,i(-i/6j-ai)+i) if s ^ o, K ^ J 

by Corollary IA.lt P ar t (b) is proved. 

(c) Let r and q% be such that the conditions of Lemma IA.1I are satisfied for all j where 
bj > —1. For all j where bj < — 1 assume that r also satisfies 

bj + r<-l. (74) 

Let Q2 and Q3 be defined as in Lemma [A. 71 for j — 1. By the same arguments employed 
in that lemma, for sufficiently small 5, if (w,z) G B$ \ (Q 2 U Q 3 ) then either ^(iu, z) > 5 
or #1(10, 2) G Q(r, 5). In the latter case, if all bj > —1 then gi t k(w, z) > 5 (for small enough 
5) for some I and k due to Corollary IA.l( b). If there exist bj < —1 then gj t \{w,z) > 5 
for small 5 due to (1741) . Therefore, gi t k(w,z) > 5 for some / and for all (117,2), such that 

We now consider (w, z) £ B$\ (Q2 U Q3). There are two (generically) mutually exclusive 
cases: 
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• Suppose that at least one of the inequalities, rii<i<m a j > 1 or a t — &t < for some t, 
is satisfied. Denote = Q 2 ^ and define the sets Q\ I' 1 and yo in the same way as 
in the proof of the part (b). Since lim^oo /ii-fcm,i(z/o) = oo, oi_ is arbitrary large. 

• Suppose that Yli<j< m a j < 1 an d a t — h > for all t. There exist q 2 > such that 
Y[i<j< m ( a j + 92) < 1- By Corollaries IA.2I and IA.3I there exist limits as k — > 00 of fi 
bounding Q\ ~ km+ ^' 1 ( e ) f or SO me e < e. Hence, finite values of 6j can be found such that 
Lemmas IA.3I and IA.4I hold true for Q^^ mJr ^ 1 for any k > 0. For any 5 < min^e^), by 
Lemmas IA.3I and IA.4I the following is satisfied 

\ 9j (w, z)\ > {w, zp +q2 for all {w, z) e Qi; 3 km '\S). (75) 

Take any (wa,za), |(iwo, ^o)| — oc < 8. There exists k > such that a kai ' m > 5, where 
Qij = Tli<s<j( a s + 92)- Hence, due to (J75]), if \gi, s (w , z )\ < 5 for all 1 < / < m and 
s < k, and \g k {w , z )\ < 5 then g k {w , zq) G B$ \ (Q 2 U Q 3 ). 

QED 



B Types B and C cycles 

In this section we present proofs of Theorems and Lemmas employed for calculation of sta- 
bility indices of Types B and C cycles. To leading order, the maps gj : — > 
associated with the cycles of Types B and C reduce to gj(w,z) = (Ew aj , Fw h3 z) and 
gj(w,z) = {Ew b: > z, Fw a i), respectively. As noted in Section I4T21 it suffices to consider only 
positive values of w and z. In the coordinates (C, r/), ( = In z and r\ = In w, the maps gj take 
the form: 

(In what follows, the constants E and F are ignored, see discussion in Section 14.21 ) The 
transition matrices of the maps are 

= (J; I) 

for cycles of Types B and C, respectively. Recall that the coefficients aj and bj of the map 
gj are related to the eigenvalues of linearisation of near £j dbS CLj — Cj /ej and bj = —tj/ej. 
As in Appendix |Aj the stability indices are calculated in terms of exponents of the maps gj, 
aj and bj. For the map g = g m o . . . og 1 the transition matrix is M = M(g) = M m ■ ■ ■ Mi. We 
introduce the notation: Mj t k and M^ denote transition matrices for the maps and g^\ 
respectively; M^' — Mj • • • Mj\ A{, A2, v{ = (v{i,v{ 2 ) and v 2 = (v ii,v J 22 ) denote eigenvalues 
and associated eigenvectors of the matrix M^\ respectively. If the eigenvalues are real, 
X{ > \ 2 is assumed. 
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B.l The set C/-°°(M) 

A necessary condition for (w,z) to belong to £>f (see Subsection 13 .31) is that g k (w,z) is 
bounded for all k by a small 5 > 0. Since g k (w,z) is bounded, in the new coordinates 
(C,T)) = (lniu, lnz) the iterates (Ck,Vk) = 9 k (C,v) are bounded from above: (k < S and 
i]k < S for some large in absolute value negative 5. Due to linearity of g, this generically 
implies that lim^oo g k ((, rj) = (— oo, — oo). 
We denote 

U-°°(M) = {(x,y): x < 0, y < 0, lim M n (x, yf = (-oo, -oo)*}. 

n— >oo 

Lemma B.l The dependence of U~°°(M) on eigenvalues and eigenvectors is as follows: 
(i) If the Xi are complex, then U~°°(M) = 0. 

(ii) If the Xi are real and X\ < 1 or \X 2 \ > X\, then U~°°(M) = 0/ 
(Hi) If the Xi are real and VnVi 2 < 0, then U~°°(M) = 0; 

(iv) If the Xi are real, X\ > 1, vnVi 2 > and v 2 \v 22 < ; then 

U-°°(M) = {(x,y): x < 0, y < 0}; 

(v) If the Xi are real, X\ > 1, Viiv u > (whereby we assume v\\ > and v u > 0), and 
V21V22 > 0, then 

U-°°(M) = {(x,y): x<0, y<0, (v n v 22 - v l2 v 21 )-\v 22 x - v 2l y) < 0}; 

(vi) If the Xi are real, X\ > 1, V\\V\ 2 = and X 2 < 1, then U~°°(M) = 0; 

(vii) If the Xi are real, X\ > 1, vuV\ 2 = 0, X 2 > 1 and v 2 \v 22 < 0, then U~°°(M) = {(x,y) : 
x<0,y<0}; 

(viii) If the Xi are real, Ai > 1, v n vi 2 = (whereby we assume Vn > and v\ 2 > 0), X 2 > 1 
and v 2 \V 22 > 0, then 

U~°°{M) = {(x, y) : x < 0, y < 0, {v u v 22 - v^n)' 1 {v 22 x - v 2l y) < 0}. 

Proof: Let the eigenvalues be complex conjugate, Ai^ = se ±l( ^ . Then 

M n (avi + /3v 2 ) = s n ( vi (a cos n<p — (3 sin ncfi) + v 2 (a sin ncj) + (3 cos n<p) ) = 

s n (cosn(j)(avn + f3v 2 i) + sinn</>(— $v\\ + av 2 i), cos n<p(av 12 + j5v 22 ) + sinn0(— /3t> 12 + av 22 )). 

Because the eigenvalues are complex, <p 7^ kir. Hence for any iVo > there exists N > Nq 
such that 

x n = s n (cosn0(afn + f3v 2 i) + smn(f>(—f3vn + av 2 i)) > 0. 
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This proves part (i). 

If the eigenvalues are real and distinct, the map M n in the basis comprised of the eigen- 
vectors Vx and v 2 , (x, y) = hiWi + a 2 v 2 , takes the form 

«•(*>,*») » i % )=^(:;:)+^(:::)- <™ 

If < 1 and | A2 1 < 1 (recall that A2 < Ai), then for any (hi, h 2 ), (!76|) has a finite limit as 
n — y 00. If |A 2 | > Ai, then A2 < and hence in ( J76l) the sign of h™ (j = 1, 2) (176]) alternates 
for odd and even n, if n is large enough. Part (ii) is proved. 

Assume that Ai > 1 and | A2 1 < X\. If t>iit>i2 7^ 0, to leading order hi = XihiVn and 
= X±hiVi2 for n — > 00. Thus, if the signs of v\\ and V12 are different, then the limits of 
hi and h^ have different signs, and (iii) is proved. If VnV 12 > (and assuming without 
any loss of generality that they are positive), the limits of h^\ j = 1,2, are —00 in the 
points (x, y) such that hi < 0. hi is negative for any x < and y < 0, if f 21^22 < 0, which 
proves part (iv). If v 21^22 > 0, then the set of (x, y) for which < satisfies the inequality 
(^11^22 — 'Wi2fii)~ 1 (^22^ — v 2\y) < 0, and so part (v) is proved. 

Assume that V12 = 0; if A2 < 1, in (1761) the limit of h 2 for n — 00 is either +00 or it 



5 2 -— , VM-JJJ/ H«G limit Ul l(, 2 

does not exist. Thus, part (vi) is proved. The proofs of statements (vii) and (viii) are similar 
to the proofs of (iv) and (v) and are omitted. QED 

Lemma IB .11 can be used to give the dependence of U~°°(M) on the matrix entries, 
M= (ay), i,j = 1,2. 

Lemma B.2 Let Ai and A2 (Xi > A2, if they are real; generically X\ 7^ A2 ) be the eigenvalues 
of the matrix M = (ay), an > «22, and vi and v 2 be the associated eigenvectors. Then 

(a) (i) the eigenvalues are real if and only if 



(an — O22) 2 . ^ n /™x 

- + ai 2 a 2 i > (77) 



4 



(ii) Ai > 1 if and only if 



1 Oil + a 22 . . > , 

max , an + ^22 - ^110-22 + ^12^21 > 1 (7t 



(iii) Xi > IA2I if and only if 

(iv) v 11^12 >0 «/ and on/?/ i/ 



> (79) 



a 2 i > 0. (80) 
(b) If Xj, j = 1,2, are real, then ai2t>22 = (A 2 — aii)^! and X2 — an < 0. 
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Proof: The eigenvalues of the matrix M can be expressed as 

Al ' 2 = 2 ± V V 2 J +ai2fl21 - 

Statements (a)(i)-(iv) follow from examination of this formula and on noting that the eigen- 
vector vx satisfies a 21 vu + a 22 v \ 2 = \vi 2 . If the eigenvalues are real, then A 2 — an < and 
v 2 satisfies a\\v 2 i + a 12 t> 22 = A 2 fn, which proves statement (b). QED 



Corollary B.l Assume that entries of a matrix M satisfy the conditions (i)-(iv) of Lemma 
M Then U-°°(M) = {(x,y) : x < 0, y < 0, (A 2 - a u )x - a 12 y > 0} for a 12 < 0, and 
U-°°(M) = {(x,y): x < 0, y < 0} /or a i2 > 0. 



B.2 Maps and neighbourhoods 

The condition w 2 + z 2 < e in (C,^) coordinates is equivalent to the condition^ max(C, i]) < R, 
where R < 0; small e corresponds to large \R\. In this subsection we examine how R- 
neigbourhoods of (—00, —00) are transformed by linear maps. 
Let M be an invertible linear map M : IR 2 — > M 2 . Define 

Ur = {(C,v) I max(C,77) < R} 

and 

U R (a 1 ,/3 1 ,qi;a2,/32, q 2 ) = 
{(C, rj) e % : ("i + ?i)C + (A + Si)?? < 0, (a 2 + q 2 )C + {$2 + q 2 )v < 0}, 

where R < 0. 

Lemma B.3 For any S < and g > there exists R < saca taai 

M(Ur n M- 1 C/ (1, 0, -g; 0, 1, -g)) C U s . (81) 

Proof: We split the neighbourhood in two parts 

U (l,0,-q;0,l,-q) = V 1 UV 2 , 

where V\ = U Q {\, 0, — g; 0, 1, — g) fl U s and V 2 = U (l, 0, — g; 0, 1, — g) \ Us and denote 

R = mm{(,r ] : ((,!))£ M"V 2 }. 

i? is finite, since V 2 is bounded and M invertible. The inclusion f lHTj) takes place for any 
R < R, because 

M(U R n M- 1 C/ (l, 0, -g; 0, 1, -g)) n K 2 = 
due to R < R. QED 

Note, that if ((,r)) G U R \ M _1 C/ , then max(C,r/) > 0, where (C, 77) = M((,r]). 

8 In other words, instead of |(u>, z)| = (w 2 + z 2 ) 1 / 2 , an equivalent norm \(w, z)\ — max(|w|, \z\) can be 
employed. 



50 



Lemma B. 4 Denote U = U (ai, ft, — q\, a 2 , ft, — q 2 ). Suppose U C U °°(M) ^ and 
{vu,v 12 ) e [/. 

(%) 7/A 2 > 0, inen M fc (f7) C U for any k > 0; 

(ii) M 2k (U) C £7 for any k > 0. 

Proof: (i) If (x,y) G U, then a(x,y) G Z7 for any a > 0. If (x,y) G Z7 is represented 
as (x,y) = avi + ftv 2 , then, due to convexity of U, arvi + /3v 2 G C/ for any \(3\, such that 
|ft < \0\ and > 0. 

For any (x, y) = av! + /3v 2 G £7 the identity 

M k (x, y) = X k (a Vl + /^/Afo) (82) 

holds. Since Ai > A 2 > 0, due to the arguments above, M k (x,y) G U. 

(ii) If is even, the identity (1521 implies the statement for negative A 2 as well. QED 

Lemma B.5 Consider the set 

ZJ-°°(M) = [/ (ai,ft,0;a 2 ,ft,0). 
For any q\ > 0, q 2 > and S < there exists an R < such that 

(i) If ^2 > 0, then M k U R (a 1 , ft, -q x \ a 2 , ft, -92) C /or any fc > 0; 
(iij M 2k U R (a u ft, -gi; a 2 , ft, -q 2 ) C C/ s /or any k > 0. 

Proof: (i) Consider a set comprised of two line segments, one segment being x — — 1, y — 
[— 1, 0] and the other one y = —1, x = [—1, 0]. Denote W = W D £/o(«i, ft, — <?i; «2, ft, —92)- 
The iterates (xk,Vk) — M k ( y x ,yo) t are bounded from above for any (x ,yo) £ W by some 
Q < 0, because W C U~°°(M). Due to compactness of W, the bounds for the iterates are 
uniform for all (x , y ) G W by some Q < 0. Since the maps M fc are linear, max(a; fc , y k ) < Q 
for a given (xo,yo) implies (xk,yk) < a< 5 f° r (^0,^0) — &(xo,yo)- Part (i) is thus proved for 
R = -S/Q. 

Statement (ii) is a consequence of Lemma IRH (ii). QED 
B.3 Main theorems 

In this subsection we prove two theorems on stability indices for maps related to heteroclinic 
cycles of Types B or C, employing the Lemmas proved in Sections IB. II and IB. 21 

Theorem B.l (reproduces Theorem \4 . 2j) Let g be a map related to simple heteroclinic cycle 
of Types B or C and Mj, 1 < j < m, its transition matrices. Suppose that for all j, 
1 < j < m, all entries of the matrices are non-negative. Then: 
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(a) If the transition matrix M = M m ■ ■ ■ Mi satisfies condition (a)(ii) of Lemma \B.S\ r then 

= oo and cr,- _ = for all j and therefore the cycle is asymptotically stable. 

(b) Otherwise, 0^4- = and Oj - = 00 for all j and the cycle is not an attractor. 

Proof: (a) Suppose that the matrix M = satisfies condition (a) (ii) of Lemma IB.2I For 
a map 

M U) : H W H i™) ) the con dition can be expressed as 

max (tr (M {j) ) , 2tr (M {j) ) - 2det (M u) )) > 2. (83) 

Hence if the condition is satisfied by M^' for any one value of j, it is satisfied for all 
1 < j < m. Any have non-negative entries, as it is a product of matrices with non- 
negative entries. Therefore, satisfies conditions (i) and (iii) of part (a) the Lemma. 
Due to the assumptions Ai > A 2 and an > a 2 2, the condition (iv) is satisfied. Hence 
U-°°(MW) = R 2 _ for all j. 

Consider the images of the lines a(— 1 + q, —q) and f3(—q, —1 + q), where < q < 1, 
a G 1R + and G R + , under the mappings M- 7 ' 1 or M J,1 M. Since all entries of the matrices 
JVP' 1 and Mi ,l M are non-negative, the images take the form a(— ri, — r 2 ) for r« > and 
a G M+. Thus, for any positive q < 1 there exists a Oj > 0, such that 

M^U (1, 0, -g; 0, 1, -q) C f/ (l, 0, 0, 1, -^) 

and 

M^MU {1, 0, -g; 0, 1, -q) C C/ (l, 0, 0, 1, -g,). 

Apply Lemma[K3]to mappings M- 3,1 and M J,1 M, setting any S < and q = qj- According 
to the Lemma, we can find Sj, such that 

M j ' 1 Us j (1, 0, -g; 0, 1, -g) C £/ s and M^MU Sj (1, 0, -g; 0, 1, -g) C (84) 

Denote 

S 1 = min Sj. 

j 

By Lemmas IB.4I and IB. 51 (where Lemma [B.5I is applied for S = S), there exists R such 
that 

M 2k U R (l, 0, -g; 0, 1, -g) C U§(1, 0, -g; 0, 1, -g) for all k > 0. 
Thus, fl84l) implies 

Af iffc E/fl(l, 0, -g; 0, 1, -g) C U s for all fe > 0. 

Hence, 01,+ > 1/g — 1 for any g, which implies = 00 and cri^ = 0. The proof holds true 
for j > 1 as well, and therefore part (a) is proved. 

For part (b), if the matrix M = does not satisfy the condition (ii) of the Lemma, 
by Lemma [Bj] the set U~°°(M) is empty, oi i+ = and <Jx - = 00. Since condition fl83l) is 
satisfied or not satisfied by all simultaneously, o"j >+ = and Uj- = 00 for all 1 > j > m. 
QED 
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Theorem B.2 (reproduces Theorem \4-3\ ). Let X be a simple heteroclinic cycle of Types B 
or C and Mj, 1 < j < m the associated transition matrices. We denote by j = ji, . . . Jl the 
indices, for which some of the entries of Mj are negative; they are all non-negative for all 
remaining j . 

(a) If at least for one of j = ji + 1 the matrix M^) does not satisfy conditions (i)-(iv) of 
Lemma \B. H then the cycle is repelling and Oj = — oo for all j. 

(b) If the matrices satisfy conditions (i)-(iv) of Lemma \B.2\ for all j — ji + 1, then 
there exist numbers (a{, P{, a\, P 3 2 ), 1 < j < m, such that 

(i) C/oW, ft, 0; 4, Pi 0) ^ 0, 1 < j < m. 

(ii) For any S < and q > there exists R < such that 

M {l ' j) {M ij) ) k {U R {a{, P{, -q; a{, Pi -q)) C U s for all I, 1 < I < m, k > 0. 

(Hi) 

\im(M^\M^) k ((, V )) = (-00,-00), for all ((, V ) E U (a{, P{ , 0; oP 2 , ft, 0). 

fc— >oo 

(iv) 

U (a{,Pi,0 ] aiPi0) = U- oo (M^)n( f| (M^)" 1 ^ n ( f| (M^ +m '^) _1 [/ 
^ If \ 2 >0 then 

U (a{,P{,0 ] aiPl,0) = U- o °(M^)n[ f| (M^)" 1 ^ 



V K/<L / \1<KL 



\1<1<L 



The cycle is a Milnor attractor. 

Proof: For (a), as noted in the proof of Theorem I4.2[ the matrices M«! will simultaneously 
satisfy, or not satisfy, conditions (i)-(iii) of Lemma IB. 21 for all j. Suppose the condition (iv) 
is not satisfied for some j = J. For any s, the iterates (M^ J )) fc M^ J,s )(x, y) on increasing k 
become aligned with (v^v^), see (ITEj) . Since v^v^ < 0, the iterates escape from Us C Hj^ 
for any S < for a sufficiently large Hence part (a) is proved. 

For (b) suppose that, for ji + 1 the matrix M^ i+1 ) satisfies conditions (iv) of Lemma [B.2 1 
The matrices Mj, ji + 1 < j < — 1, have positive entries, v- 7 = Mj_i . . . Mj l+2 Mj l+ iv^ +l , 
therefore v il v > implies f n > for any j, ji + 2 < j < Hence, it suffices 
to check condition (iv) for j = ji + 1, 1 < I < L. 

Denote 

t/,- = E7-°°(mW) n ( p| (m^'^Uo j n ( p| (M^ +m ^)- 1 f/ 

\KKL / \1<«L 
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The set is non-empty, because it includes a neighbourhood of the point (v 3 n ,v J 12 ) on the 
plane (since this point belongs to all sets in the intersection). Since all the sets are of the 
type U (ai, (3 1 , 0; a 2 , (3 2 , 0), the intersection is also of the required type U (a{, (3{, 0; a 2 , (3 2 , 0). 
Consider U = U\. 

Due to Lemma [B. 41 and definition of the set U, 

M j>k U c tr°°(M w ) for all 1 < j < m, k > 0. 

By the same arguments as employed in the proof of Theorem 14.21 this inclusion implies that 
for any q > and S < there exists R < 0, such that 

M jtk U R (a{, ft, -q; a\, ft, -q) C U s for all 1 < j < m, k > 0, 

and therefore — oo < o\. The proof for Oj with j > 1 is similar. Finally, by Theorem 12. 3[ X 
is a Milnor attractor, since the inequality — oo < Oj is satisfied for all j. QED 
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